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MEASURING STREAMFLOW UNDER ICE CONDITIONS®,1 


| A. M. Moore,” A.M., ASCE 
| (Proc. Paper 1162) 


| ABSTRACT 


| Effects of ice formation on stage-discharge relation are explained briefly. 
Decurrence and effect of surface, frazil, and anchor ice is described and re- 
corder graphs registering these effects are shown. One investigation of “si- 
dhon action” is detailed. Factors affecting the accuracy of ice-affected rec- 
ords are listed and general appraisal made of accuracy of such records. 


INTRODUCTION 
: 
i Engineering literature contains few references concerning the effect of ice 
on streamflow measurements. In 1922, J. C. Stevens presented his paper 
‘Winter Overflow from Ice Gorging in Shallow Streams”,® and in 1931 the 
society published the final report of the Committee of Power Division under 
he title “Ice as Affecting Power Plants”,4 At the Western Interstate Snow 
survey Conference in Seattle in 1940, Walter J. Parsons, Jr., presented a pa- 
ser on “Ice in the Northern Streams of the United States”.°,6 Those papers 
lescribe or discuss the formation of ice in streams but do not discuss how it 
iffects the measurement of streamflow. 

The water-resources investigations of the Geological Survey necessarily 
nclude measurements of streamflow throughout the United States under all 


Note: Discussion open until July 1, 1957. Paper 1162 is part of the copyrighted Jour- 
nal of the Hydraulics Division of the American Society of Civil Engineers, Vol. 83, 
‘No. HY 1, February, 1957. 
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conditions of weather. Thus it was logical for the Survey to study the effect — 
of ice in open channels. Water-Supply Paper 337, “The Effects of Ice on : 


Stream Flow”, by W. G. Hoyt in 1913, describes methods and techniques de-_ 
veloped up to that time for determining stream discharge under ice condi- 
tions. Since then there has, of course, been advancement in our knowledge of 
behaviour of streams during périods of ice effect. This paper summarizes 

current practices of the Survey, several manuals prepared for use within the 


Survey, and personal experience. 


Statement of Problem 


The problem of determining streamflow under ice conditions is to achieve 
a procedure for determining the daily discharge even though the stage-dis- 
charge relation is ineffective owing to ice at the control section. 

For those who may be unfamiliar with the principles involved it may be 
explained that, when ice on the stream is not a factor, most stream-flow rec- 
ords are derived from a simple correlation between stage and discharge at 
each gaging station. A water-stage recorder provides a continuous graph of 
stage, which is converted to discharge by means of a stage-discharge relatior 
known as a rating curve. Such a curve is developed for each station from 
discharge measurements made by current meter at various rates of flow. 

So far as is practicable, gaging stations are established at points where 
the stage-discharge relation will remain stable. The first riffle, or steepen- 
ing of water-surface profile, downstream from the gaging station is known as 
the control section. At this point the cross-sectional area of the stream is 
constricted; causing a tendency for the water upstream from the control to be 
ponded. Therefore, the control is the determining factor in the stage-dis- 
charge relation. Boulders, gravel, or sand washing onto or from the control 
may cause changes in the stage-discharge relation. Ice forming on the con- 
trol may cause further ponding, or backwater. 

At some stations, particularly on large streams, the controlling feature is 
the bed and banks of the river for a considerable distance downstream from 
the gage. That is known as channel control. Here, too, the formation of ice 
can cause backwater. Ice may clog the channel, or surface ice may increase 
the frictional resistance to flow with a decrease in velocity that results in a 
compensating increase in stage. Sometimes when there is a combination of 
channel control and section control, surface ice can cause backwater even 
though ice has not formed on the section control. 

Although the control may be permanent under open water conditions, ice 
forming at the control causes a temporary change in the stage-discharge re- 
lation usually by raising the gage height needed to pass a given quantity of 
water. Figure 1 compares a stage-discharge relation defined by open-water 
current-meter measurements with the plotting of measurements affected by 
backwater from ice. 


Types of Ice Formation in Streams 


Surface Ice 


of the three types of ice that form in streams, the most familiar and ob- 
vious is surface ice. As the name implies, it forms on the surface, either a: 
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Fig. 1. Rating curve for Menominee River near Pembine, Wis. 


fringe of “shore ice” or as a continuous cover spanning the stream from 
ank to bank. Surface ice that completely covers a small stream may be 
wrmly anchored at either shore. Then the amount of backwater can vary 
reatly with stage—for example, there may be no backwater effect if the stage 
lls and the ice forms a bridge above the water. On large streams, surface 
'@ usually floats so that the backwater effect is relatively constant with 

‘age but changes as the ice thickens or thins. 

Figure 2 shows the water-stage-recorder chart for a gaging station as a 
| Bnation of surface ice begins to cause a backwater effect. In this example, 
aily mean discharge remained about the same as before the freeze-up, but 
adoubtedly the momentary discharge fluctuated somewhat during each day. 
| Surface ice can cause “siphon action” also. This effect, which is opposite 
' that of backwater, will be covered more fully later in this paper. 

i] 
: 
‘razil Ice 


| “Frazil” is a French-Canadian term for finely spicular ice particles 

thich, owing to wind or turbulence, do not unite to form surface ice and 

ich, therefore, float downstream. Commonly the spicules or needles of 
group together into a slushy mass which moves downstream until it lodges 
cold rocks or shores or freezes to the underside of surface ice, causing 
ekwater, These floating masses also frequently interfere with the operation 
/acurrent meter. 


chor Ice 


The third type of ice formation, “anchor ice,” usually is an accumulation 

| spongy ice on the stream bed. When this spongy ice forms at a control 
ction it will produce characteristic “humps” or “anchor-ice rises” in the 
ige-height graph. Only rarely does anchor ice form as a solid ice cap on 

| merged rocks, and when it does the backwater effect is not large. Anchor 
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TYPICAL RISE AS COMPLETE ICE COVER FORMS 
NORTH FORK CROOKED RIVER BELOW DEEP CREEK, OREG. 
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Fig. 2. Typical rise as complete ice cover forms. 


ice results from loss of heat by radiation from the streambed to the atmos- 
phere and therefore does not form when the stream is completely covered 
with surface ice. Generally it forms during the late evening of cold, clear 
nights and dissipates when the sun’s rays warm the rocks the next day. The 
usual duration period for anchor ice is 8 to 15 hours, but if the weather is 
cloudy and cold the anchor ice may remain in place for a considerably longer 
time. Figure 3 shows a typical effect of anchor ice on a water-stage record- 
er record. The distinguishing feature of the “anchor-ice hump” is that the — 
rise is slow compared to the fall, whereas an actual increase in streamflow 
would occur in the opposite sequence, or at least the rise would be as rapid ~ 
as the fall. Therefore, simple diurnal anchor-ice rise is readily recognized 
and the graph can be corrected simply by “cutting out” the hump in the 
record. ; 
On figure 3, the dashed line represents the estimated gage height corre- 
sponding to actual discharge. The dashed vertical lines indicate noon and the 
full vertical line, midnight. Note that the discharge and the effective gage 
height tend to diminish as the backwater increases, and then increase sharp- 
ly as the anchor ice leaves the control and the detained water is released. 
Note, also, the small rises in the late afternoon. These probably are the re- 
sult of increases in discharge from releases of anchor ice upstream com- — 
bined with some melting of snow and ice during the warmer part of the day. 
Anchor-ice dams at successive riffles on a stream can detain a relatively 
large volume of water and cause sharp drops in streamflow for short period: 
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| TYPICAL ANCHOR ICE RISES 
NORTH FORK BEAVER GREEK NEAR PAULINA, OREG. 
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Fig. 3.. Typical anchor ice rises. 


: 
: 
' 
| Siphon Action 

| As has been mentioned, surface ice can cause an effect opposite to that of 
jackwater—that is, the flow is more than would be expected for the observed 
age height. This phenomenon of “siphon action” occurs infrequently, at only 
few stations, and only when complete ice cover extends for some distance 
bove and below the control section. It occurs when the hydraulic environ- 
nent resembles that of a siphon spillway in that the ice cover prevents the 
eady entrance of air. As water siphons over the control, the stage may drop 
rom that of an antecedent backwater condition to a level below that corre- 
ponding to discharge under an ice-free condition. Once the stage drops suf- 
ciently so that air can enter, the siphon is broken and the stage rises until 
ir is again prevented from entering under the ice cover and the siphon is 
gain primed. The duration of each period of siphon action may range from 
bout a minute to more than an hour, depending on the thickness of ice cover, 
mount of flow, and size of pool above the control section. Siphon action is 
ost common on artificial controls or weirs but can occur on some natural 
ontrols. 

The author and J. V. Bagley of the Boston office of the Geological Survey 
laborated on a study of siphon action at Oyster River near Durhan, N. H., 
1 January 10, 1947. During that day the temperature ranged from - 190F. to 
17°F, Briefly, the highlights of that study are as follows: 

The control (Fig. 4) was a V-shaped concrete weir about 26 feet wide situ- 
ed 30 feet downstream from the gage. Figure 5 shows this control section 
irtly ice covered, and figure 6 shows it about as it was on the day of the 
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Fig. 4. Control for Oyster River near Durham, N. H., under open- 
water conditions. 


Fig. 5. Control for Oyster River near Durham, N. H., partly 
frozen over. 


study—that is, completely covered with ice and with siphon action occurring. 
Figure 7 shows the recorder graph for the study period. The slight fluctua- 
tions in the graph beginning about 11:30 a.m. on January 10 were caused by 
engineers clearing a measuring section about 80 feet downstream. At this 
time the stage was actually slightly below the point of zero flow for open- 
water conditions. A discharge measurement was made which showed the 
flow to be 7.02 cfs. The siphon was then broken by chopping a hole in the ice 


Fig. 6. Control for Oyster River near Durham, N. H., completely 
frozen over. 
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Fig. 7. Effect of siphon action. 


just upstream from the weir and the stage rose to a gage height of 0.24 foot 
before the hole froze over. As siphon action began again the stage dropped 
to about 0.02 foot, at which point all ice was chopped from the control and the 
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stage rose and leveled off at a gage height of 0.39. The rating curve, which 
was well defined, showed that the discharge for that stage was 5. 4cfs. Peri- 
odically throughout these tests the velocity of the water at a fixed point in the 
cleared section was measured by current meter. The velocity was 0.99 fps 
before the siphon was broken, dropped as low as 0.30 fps after the siphon was 
broken, and then increased to 0.56 fps. When the hole froze over and siphon 
action began again the velocity jumped to 1.50 fps and then decreased to 0.83 | 
fps just before the control was completely cleared. 
From that study it was concluded that: 


1. As “make and break” siphon action continues, velocity and discharge 
are continually changing, rising to a maximum as siphon action begins and 
dropping to a minimum as the siphon is broken (could be zero if the pool had 
been drawn below the zero-flow point). 


2. A single discharge measurement is of little or no help in determining 
mean discharge unless it is made after the siphon is broken and flow has 
reached equilibrium, or is made above the gage pool. 


Factors Affecting the Accuracy of Streamflow 
Records Under Ice Conditions 


Frequency and Thoroughness of Gaging-Station Inspections 


Nothing can take the place of on-the-spot determinations of actual flow and 
ice conditions. During open-water periods, inspections usually are made 
once every 4 or 5 weeks, but during the winter, inspections are made twice a 
month at some stations. At regulated stations, such as those below power- 
plants that carry a variable load, it may be necessary to make two measure- 
ments during each winter visit—one at the high stage of the regulation and 
another at the low stage. The backwater effects from ice may be markedly 
different at these two stages. 

Measurements may vary widely in accuracy owing to adverse ice condi- 
tions—for example, ice floating and jamming in the stream, or even ice and 
water in alternating layers. Also, extremely cold temperature may cause 
meter trouble—for example, ice may form on the meter as it is lifted from 
the water to move it from one hole in the ice to another. Sometimes it is 
necessary to have at hand a pail of heated water in which the meter can be 
immersed while the engineer moves from one hole in the ice to another. Re- 
gardless of the precautions taken, some measurements are unavoidably less 
accurate than those made at the same site under ice-free conditions and at 
moderate temperature. 

Accuracy of adjustments for ice effect between inspections are greatly im- 
proved if the inspecting engineer makes complete and thorough notes about 
ice conditions. Not only does his measurement give the flow and the amount 
of backwater at the time of his inspection, but his notes give dependable clues. 
to the amount of backwater effect at lower or higher stages. Some engineers 


make a practice of photographing the control section and its ice cover during: 
each inspection. 


Type of Ice Formation 


Surface ice can cause much uncertainty as to the discharge because its 
effect does not show clearly on the recorder chart. For example, it may be 
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evident that a backwater effect exists and even that it is increasing or de- 

| creasing, but the amount of that effect cannot be determined directly from the 
| recorder chart. Anchor-ice rises, however, are clearly recognizable on the 

| recorder chart, and their effect usually can be corrected directly on the chart 
} as shown in figure 3. The shape and position of the rise just before noon is 

, estimated on basis of similar rises at gaging stations not affected by back- 

_ water from ice. Some loss in accuracy occurs in all records affected by sur- 
_ face ice—particularly if siphon action occurs, but records affected by anchor 

) ice alone are nearly as accurate as the open-water records. 
| Type of Gage 


A water-stage recorder provides a much better winter record than a staff 
| or wire-weight gage, because the recorder graph itself usually provides de- 

/ pendable clues to the discharge or at least to discharge trends. For several 
years the Geological Survey has operated some thermograph recorders that 

give a record of stream temperature, usually within 19°F. If the temperature 
| Fecord shows that water is more than 1° above freezing it is safe to assume 
, that ice is not forming, although ice that had formed previously could still be 
| Causing backwater. The water-temperature record is especially useful in 

| fixing the time that ice begins to form. 


| Comparison with Records for other Gaging Stations 


When surface ice is present the chief basis for determining the probable 
| discharge between discharge measurements is usually a study of runoff by 
) comparison with records of other nearby streams. Even though the record 
used for comparison may be ice-affected also, it affords a different and inde- 
| pendent set of base data—another recorder chart and a different set of cur- 
_rent-meter measurements. Without a nearby record that compares well with 
_the record in question, the accuracy of the computed discharge may be great- 
ly reduced. Fortunately the network of stations is sufficiently extensive that 
ordinarily more than one record can be used. However, this is not always 
i true, especially in the more sparsely settled areas. For example, in south- 
eastern Oregon, stations are quite widely scattered and, owing to diverse hy- 
drologic environment, the records do not compare well with one another. 
Comparisons always should be used with caution, as the relationship be- 
tween the flow of two streams may vary considerably during the year. Figure 
8 shows the usual method of plotting discharge records with the time scale 
horizontal and with the discharge scale vertical and logarithmic. The solid 
line represents the apparent discharge obtained by applying the recorded gage 
height to the stage-discharge relation curve and the cross on January 26 rep- 
resents the apparent discharge at the time of the measurement. The dashed 
line represents discharge after the effect of backwater was removed. Full 
basis for the adjusted.graph is not shown here. Figure 9 shows the tempera- 
ture graph and independent discharge record that also formed part of the basis 
for the adjusted record. 

Useful comparison does not require that the shapes of the hydrographs be 
exactly alike as long as the trends are similar. On figure 9, note that the two 
streamflow records and the concurrent weather record respond to the effect 
of temperature during the winter period. The record for John Day River at : 
Service Creek, Oreg., is unaffected by backwater from ice, and for eee y; 
the adjusted record only is shown for North Fork John Day River. In actu 


1162-10 


HY 1 


NORTH FORK OF THE JOHN DAY RIVER AT 
MONUMENT, OREG. 


4,000 


| 


ENGINEER'S 
NOTES 


=: 


O= ane EASURE 


a 


A, CONTROL 
CLEAR 


B. BANK ICE, 


SLUSH ICE, 


+ 


NO EFFECT 


WATER 33°F. 


C. BACKWATER 


FROM ICE 
COVER 


DOWNSTREAM 
WATER 32 °F. 


100 


Fig. 8. Adjustment for backwater from ice. 
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Fig. 9. Comparison of daily discharge. 


practice all records are usually plotted on separate hydrographs, and, to 
make a comparison, one is placed over the other on a light table. At times a 
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streamflow record can be compared effectively with a powerplant record ex- 
pressed in units of either discharge or power output. As a rule, powerplant 
operators maintain fairly accurate records, but even those records that may 
not be completely accurate quantitatively are useful in showing trends. 


Comparison with Weather Records 


Good records of air temperature and precipitation are an invaluable aid in 


| making adjustments for ice effect. The temperature record helps the engi- 


neer to decide whether ice is forming, increasing, or dissipating, and whether 
precipitation is rain or snow. 

For stations which do not have a nearby discharge record for comparison, 
and for which the recorder chart does not provide a dependable clue, it may 
be necessary to adjust for ice effect almost solely on the basis of weather 
records and a few measurements of discharge. Discharge usually follows re- 
markably well the “ups and downs” of the aid-temperature record, and oc- 
casional discharge measurements help fix the estimated rises and drops 


| within reasonable limits. 


Comparison with Ground-Water Records 


During the winter, when much of the precipitation occurs as snow, a major 
part of streamflow is derived from ground water. A record of fluctuations in 
ground-water level for such an area can be very useful in estimating 
discharge. 


Thoroughness and Skill of Office Engineer 


Computation of discharge from a gage-height record affected by backwater 
from ice is a painstaking job that involves much detail. The experienced en- 
gineer overlooks nothing and makes the adjustments on a day-by-day basis, 
giving consideration to recorder charts, to all pertinent discharge measure- 
ments, to engineers’ field notes and photographs, and to records of nearby 
streams, ground-water levels, and weather. 


Accuracy of Ice-Affected Records 


Accuracy Statements in Water-Supply Papers 


When only diurnal anchor ice is present, the accuracy of streamflow rec- 
ords generally is not affected. When surface ice is present for many days 
the records of streamflow for those days must be rated lower in accuracy, 
but usually not more than 10 percent below the accuracy of the open-water 
records for the same station and commonly by no more than 5 percent. If 
the “Remarks” paragraph of the station description in a water-supply paper 
states “Records excellent except those for periods of ice effect which are — 
good,” the accuracy of daily discharge is considered to be, in general, within 
5 percent of open-water conditions and within 10 percent for icing condone: 
A qualification of “fair” indicates accuracy within 15 percent, and poor” in- 
dicates that the errors could be more than 15 percent. A qualification of 
“poor” may be necessary for ice periods if the only clues to the discharge 
are weather records and a few measurements spaced a month or more apart. 
With our present knowledge concerning ice formation, and with our present 
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methods for adjusting records, it seems doubtful that any ice-affected record 
should be more than 25 percent in error. Most such records should be less 
than 15 percent in error, or “fair.” 


Present-Day Records Versus Older Records 


Some surprisingly good adjustments for ice effect were made in early 
years but, in general, these older adjustments are less reliable than those 
made in recent years, as most early records were based on gage heights 
from nonrecording gages. The most significant advance in knowledge of 
stream behavior under ice conditions that has occurred since Water-Supply 
Paper 337 was written is a result of the development and more general use of 
the water-stage recorder. Also, in the early days of stream gaging, engi- 
neers did not have the time to make the detailed examinations which are 
commonplace today. At the present time, streamflow records the Nation 
over are being reviewed and compiled in summary form—specifically, all 
records through 1950 are being carefully scrutinized and, if necessary, re- 
vised. This should do much toward correcting any large errors in adjust- 
ments for ice effect during early years. 
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J. G. Slater,! J. R. Villemonte,2 A.M. ASCE, and H. J. Day3 
(Proc. Paper 1163) 


SYNOPSIS 

| This paper presents the results of an investigation conducted at the Hy- 
' draulics Laboratory of the University of Wisconsin on the energy loss in oil 
‘flow in straight, smooth pipes at line pressures up to 2000 psi and tempera- 
‘tures up to 120° F. A total of 1200 tests were completed using Klondyke light 
| > medium hydraulic oil and fuel oil No. 2 in four seamless steel pipes, 1/4, 
1/2, 1 and 2 inches in diameter. The flow was varied to produce a change in 
| Reynolds Number of 146 to 135,000 for a range of atmospheric viscosities of 

0.00003 to 0.0019 sl per ft sec. The viscosities were measured independently 
“by two methods: (a) Stoke’s Law falling--sphere, and (b) Hagen-Poiseuille 
| laminar flow pipe. Both viscometers were operated at the pressure and tem- 
| perature ranges of the tests. 
There was a small consistent increase in the energy loss due to pipe fric- 
) tion at 2000 psi over that at atmospheric pressure. In the absence of signifi- 
cant experimental errors, it is concluded that the performance character- 
istics of the Stoke’s Law and Hagen-Poiseuille Law viscometers are not 
identical because of the different operating conditions of boundary, velocity 
distribution, character of distortion rate, and body and surface forces. It is 
recommended, therefore, that the viscometers be matched with the flow prob- 
len, insofar as the operating conditions are concerned, and particularly so 
when the pressures exceed 3000 psi. 
__ The standard pipe friction theory is valid for wide ranges in pressure and 
‘temperature, provided the appropriate effects of pressure and temperature 
on the fluid properties are known. For all practical purposes, the Gunaji- 
Villemonte viscosity correlation constants can be combined with the Hagen- 
Poiseuille and Blasius formulas for petroleum oil flows at Reynolds Numbers 


less than 150,000. 
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INTRODUCTION 


Modern hydraulic control circuits and servo-mechanisms use medium vis- 
cosity oils in seamless tubing under line pressures varying from a few hun- 
dred to several thousand lbs. per sq. inch. In design work it is customary to 
calculate the head loss due to pipe friction by solution of the Darcy-Weisbach 
equation 


beet f 
hy Sites or (1). 


where f is a dimensionless friction factor dependent on the Reynolds Number 


R = Ov (2) 


and the relative roughness, k/D, in turbulent flow. In oil hydraulics it is not 
unusual for the flow regimen to be laminar, in which case the head loss is de- 
fined by the Hagen-Poiseuille Law 


work Papi Ba" (3) 
= 2 5Y 


which can be rearranged in the form of equation 1 to give the relationship 


wed (4) 
Tae 


In turbulent flow, as long as the pipe wall protuberances do not disturb the 
normal relationship between the thickness of the laminar sublayer and R, the 
relationship between f and R will follow that for the so-called “smooth” pipe. 
Most oil-hydraulics circuits are operated at R < 150,000, and use seamless 
pipes whose absolute roughnesses are substantially less than the thickness of 
the laminar sublayer. In general, therefore, these are smooth pipes which 
have friction factors independent of the relative roughness. 

For turbulent flow in smooth pipes, with R < 150,000, the Blasius formula 


f = 0.3164 R°* 6) 


is quite satisfactory for low pressure flows; and for R >150,000 the Karman- 
Prandtl formula 


oye - 
qr 7 2 L0G (RAF) - 0.8 (6) 


can be used with more confidence than ease. The Blasius formula for head 
loss can be obtained by solving equations 1 and 5 simultaneously to give 


0.25 7 1.75 
h, = a ae a) 
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The hypotheses and assumptions from which these equations have been de- 
jveloped might be considered as the present-day standard pipe friction theory. 
NThe theory presupposes that f is a function only of R and can be evaluated for 
ithe given conditions of velocity, diameter, pressure and temperature only if 
jthe effects of pressure and temperature on the fluid properties are known. 
iThe designer of oil-hydraulics systems is confronted, therefore, with two 
real questions: 


(1) Does the standard pipe friction theory apply to present-day high pres- 
sures and temperatures? 


(2) What are the effects of pressure and temperature on the properties of 
the fluid in question? 


| In regard to Question 1, a search of the literature reveals that there is 
ittle or no experimental work on friction losses in pipes at high pressures 
ithat would serve to confirm the validity of the theory, except that some in- 
conclusive tests have been made with gases at pressures up to 4,400 psi4,9,6, 
philosophical answer to this question, however, would be that it is logical 
vand reasonable to expect it to be applicable, but there seems to be no prima- 
facie evidence one way or the other. 
| In regard to Question 2, a general correlation of the effects of pressure 
wand temperature on the properties of gases is well enough defined for most 
sengineering problems, but such a correlation for liquids is not available. 
| In 1947 an investigation was initiated at the Hydraulics Laboratory of the 
ABiversity of Wisconsin to study the answers to the above questions. The 
project was necessarily divided into two parts, (a) the study of the pressure- 
PP iccasite torap irae relations of five petroleum oils, and (b) the 
study of energy loss in oil flow in straight, smooth pipes at high line pres- 
ure and a modest range in temperature. V. N. Gunaji and J. R. Villemonte? 
have completed Part (a) of the project and reported their findings, a summary 
pf which is presented in the Appendix. 
| This paper describes the work and findings in connection with Part (b) of 
the project, which was started in 1948 and completed in 1952. 
| 


Nomenclature 
The following letter symbols are adopted for use in this paper and have 
Himensions expressed by any compatible system of units. 


| 
| Aj = area of pipe orifice. 
| C =coefficient of discharge for pipe orifice. 


i. Wildhagen, M., “Uber den Stromungswiderstand hochverdichteter Luft in 


| Rohrleitungen”, Zeitschrift fuer Angewandte Mathematik and Mechanik, 
| vol. 3, 1923, page 181. 

. Newitt, D.M., and Sirkar, S.K., “The Flow of Gases at High Pressures 
through Metal Pipes”, Transactions, Institution of Chemical Engineers 
(London, England), vol. 9, 1931, pages 63-73. oe 

. Moulton, R. W., and Beuschlein, W.L., “A Study of the Flow of Air in 

Tubes in the Pressure Range of 1 to 300 Atmospheres”, Transactions, 
American Institute of Chemical Engineers, vol. 36, 1940, pages 113-133. 
Gunaji, V. N., and Villemonte, J. R., “Effect of Pressure and Temperature 
on the Viscosity of Petroleum Oils”, Proceedings of the Third Midwest- 


ern Conference on Fluid Mechanics, 1953, pages 719-740. 
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D =diameter of pipe. 

Do = diameter of pipe orifice. 
f pipe friction factor (Darcy-Weisbach). 

g = gravitational acceleration. 

h = pressure head differential for pipe orifice. 

he = head loss due to pipe friction. 

k = pipewall roughness: average height of protuberance. 
L = length of pipe. 
Q 

R 

Vv 

v 


now 


= volumetric discharge rate. 
= Reynolds Number of orifice. 
= mean velocity in pipe. 

= kinematic viscosity. 


Experimental Apparatus 


The test pipes were made from seamless high carbon steel tubing, ASTM 
Spec. A 106-48T, the dimensions of which are shown in the following table. 


Nom Sche- Actual Test Test Entr Exit Test 


Dia dule Inside L L L L D 
Dia. D D D Kk 
1/4 80 0.2988 3.C0 120 120 30 6000 
1/2 40 0.6187 7.00 136 116 29 5600 
1 40 1.049 11.00 126 103 aly y 17500 
2 80 1.952 19.99 123 117 9 13000 


The average height of pipewall protuberance, k, was measured by longi- 
tudinal profilometer readings on sections of the same material that was used 
to make the test pipes. The flanges at either end of the test section were 
welded so that the flange faces were flush with the pipe ends and the pipes 
concentrically oriented. Four radial piezometer holes of 1/8 inch diameter 
were drilled in the flange immediately preceding the joint at either end of the 
test section. These openings (8 for the 2-inch test pipe) admitted static pres 
sures to a common collector ring for pressure measurement. 

A schematic diagram of the test circuit is shown in Fig. 1. The main loo] 
was fabricated from 2-inch schedule 80 seamless steel tubing. Flow rates 
up to 100 gpm were circulated by a Northern Ordnance Company herringbone 
gear pump driven by a 10-HP variable speed motor. The bypass provided 
an additional means for controlling the discharge rate. Pressure in the line 
was initially created by a Denison hydfaulic test stand and measured by a | 
0-3000 psi Bourdon gage. The stand’was then disconnected because the Gree 
hydro-pneumatic accumulators were able to compensate for the small leak- 
age. Fluid temperatures were observed using nickel wire resistance-bridge 
thermometers immersed in the fluid. The discharge rate was measured by ; 
flange-tapped pipe orifice connected to a Statham resistance-bridge strain 
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gage with a range of 0-40 psi. Similar gages with ranges of 0-5 and 0-100 
psi were used to observe the pressure drop due to pipe friction in the test 
section. The resistance-bridge temperature and pressure apparatus was 
operated with d.c. voltages in conjunction with high-sensitivity microamme- 
ters for measuring the out-of-balance. 


Laboratory Techniques 


Figures 2 and 3 show the effect of temperature and pressure on the prop- 
erties of the two test fluids, Klondyke light-medium hydraulic oil and No. 2 
fuel oil. This information was taken from a study reported by V. N. Gunaji 
and J. R. Villemonte8, a brief summary of which is in the Appendix. The 
specific gravity of the oil at atmospheric pressure was measured by means 
of standard pycnometers. The effect of pressure was computed from the 
formula of Dow and Fink’, equation 10. Gunaji obtained the molecular vis- 
cosity of each oil by applying Stoke’s Law to the motion of a falling sphere it 
a pressure cell and accounted for the shape and boundary factors through 
comparison of his atmospheric data with calibrated Cannon-Fenske-Ostwald 
viscosity pipettes. 

The electric strain gages were calibrated at frequent intervals by dead- 
weight methods under-atmospheric line pressure. A water column was used 
for the 0-5 psi gage and the others were loaded by a standard dead-weight 
gage tester. The latter calibrations were checked and verified dynamically 
by installing the gages in the circulating system in parallel with oil-mercur: 
differential manometers for tests at low line pressures. Tests conducted by 
the Statham Laboratories have indicated that the gage calibrations are inde- 
pendent of the line pressure. 

The four orifice openings were sized so that the differential strain gage 
connected to the pipe orifice had an output always in excess of 10 microam- 
peres on the 0-100 microammeter. The orifices were calibrated in place by 
rearranging the circulating system slightly for discharge into a weighing 
tank. Tests were run with Klondyke light medium hydraulic oil and the co- 
efficient of discharge was calculated from the orifice equation 


Q =CA,NZ2gh ( 


and plotted as a function of Reynolds Number, defined as 
poe D.NZgh ( 
0-9 


The resistance thermometers were calibrated by immersion in a con- 
trolled temperature oil bath for comparison at temperatures of 54° to 120° 
F with a calibrated, precise, mercury thermometer reading to 0.1° F. 

The accumulators were precharged prior to a test series for a given line 
pressure. The test stand was used to maintain a line pressure of approxi- 


8. Ibid. 
9. Dow, R. B., and Fink, C. E., “Computation of Some Physical Properties 


Lubricating Oils at High Pressure”, I-Density, Journal of Applied Physi 
vol. 2, 1940, pages 353-357. 
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mately 200 psi during the initial bleeding to clear the lines and gages of air 
first statically, and then with the circulating pump in operation. After build 
ing up a pressure in excess of that desired for the test series, the test stand 
was disconnected, and bleeding continued until the line pressure decreased 
to that desired for the test. The strain gages were then balanced by stoppi 
the oil flow and adjusting the zeroizing circuit if necessary. This was re- — 
peated after each set of three runs, and the entire system was bled after : 
every six runs. The line pressure was recorded prior to each test. Approxi 
mately 40 seconds elapsed during the observation of the data which consisted 
of the three simultaneous sets: 


1. head loss in test section, oil temperature downstream. 


2. head loss in test section, oil temperature upstream. 
3. pressure drop in orifice, oil temperature downstream. 


At the end of each testing period, an oil specimen was drawn from the cir 
cuit and its viscosity checked. This verified the constancy of the test fluid 
composition. 


Scope of the Tests 


A total of 1200 test runs was completed using two petroleum oils (Klon- 
dyke light-medium and fuel oil No. 2) in four pipes, 1/4, 1/2, 1 and 2 inches 
in diameter. The Reynolds Number, R, ranged from 146 to 135,000 in each 
of three pressure series of 150, 1000 and 2000 psi, for an atmospheric vis- 
cosity range of 0.00003 to 0.0019 sl per ft sec. The temperature range of the 
tests was about 50° to 120° F. 


Results 


The test data were analyzed by substitution of the observations in equa- 
tions 1 and 2 in order to obtain simultaneous values of f and R. It was not 
necessary to correct for expansion of the pipe due to pressure because 
analytical studies gave a diameter increase of only 0.03 per cent under the 
worst combination of conditions. The pipe orifice calibrations were substan: 
tiated by the consistency of data in regions where there was an overlap of 
results from two different sized orifice openings. In general, there was goo 
agreement among the pipe friction tests on all sizes of pipe at all pressures 
for both oils. 

Figure 4 shows the variation of f with R on log-log paper for the data on 
all test pipes with both fluids at all pressures. The reference line in the 
laminar region is equation 4; in the turbulent region, the Blasius formula, 
equation 5. Figure 5 shows the above information for the 100-200 psi line 
pressure only, Figure 6 for 1000 psi, and Figure 7 for 2000 psi. It may be 
seen that the data are quite consistent for the Klondyke light-medium oil. _ 
There is, however, a greater scattering of the tests made with No. 2 fuel oil 
because for these tests simplified test procedures were used to expedite the 
work. A detailed study of the more consistent Klonkyke data follows. 

To facilitate comparison, straight lines were fitted to the log-log data by 
the method of averages. Figure 8 shows the ratio of the empirical equations 
for 100-200 psi data to those of Hagen-Poiseuille and Blasius. Figure 9 giv 
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the ratio of the 1000 and 2000 psi data to the 100-200 psi data. On the basis 
of the latter comparison, it is found that the 1000 psi lines check the low 
pressure data within the limits of accuracy for the tests. On the average, th 
empirical equation for laminar flow at 2000 psi is 4.6 per cent higher than 
the 100-200 psi data, and for turbulent flow the value of f is 1.8 per cent 
higher. 

The above deviations must be due either to a combination of experimen- 
tal errors or the inability of present-day pipe friction theory to account for 
all aspects of the mechanism of fluid resistance. In studying the experimen- 
tal errors, the basic measurements of discharge rate, pressure difference, 
temperature, density, and viscosity which were observed at low pressures 
must be compared with those at high pressures. Errors in mass density are 
eliminated immediately because the total density change, by itself, is less 
than one per cent. Errors in the measurement of pressure differential 
across the orifice and in the test sections are negligible because the manu- 
facturer conducted a most extensive series of pressure tests on their strain 
gages and found the pressure effect to be insignificant. The orifice was cali- 
brated at low pressures and the results used for high pressures as well be- 
cause the coefficient of discharge is a function largely of contraction rather 
than energy loss phenomena. A substantial change in the orifice energy loss, 
therefore, would produce a very small change in the discharge coefficient. 
The operation of the resistance-bridge thermometers is not affected by 
changes in pressure. 

The Newtonian viscosities used in this work were measured by Gunajil0 
using a new type high pressure Stoke’s Law falling-sphere viscometer. It 
was calibrated at atmospheric pressure by comparison of its performance 
with those of precision Cannon-Fenske-Ostwald pipettes. The calibration 
should be unaffected by pressure changes within the cell, provided appropri- 
ate corrections are applied to the sphere size and elastic and density prop- 
erties of the sphere and oil. Such corrections were made with great care anc 
are believed to be accurate. The tests show, however, that the Newtonian 
viscosities measured by the Gunaji viscometer at 2000 psi are uniformly 4.6 
per cent lower than corresponding ones computed from the laminar flow pipe 
friction data using the Hagen-Poiseuille Law, provided the other errors are 
assumed to be negligible. This is deduced from the fact that f and v have a 
direct linear relation in laminar flow. It is interesting to note that a 4.6 per 
cent change in viscosity in turbulent flow would produce a 1.2 per cent change 
in f, since f and v are related according to the Blasius formula. This shift 
in f compares quite favorably with the observed change of 1.8 per cent men- 
tioned previously. 

These deviations, therefore, seem to be caused by the different perform- 
ance characteristics of the two types of viscometers. The assumption that 
they might give slightly different results seems valid since (a) the boundary 
structures are essentially inverted, (b) the resulting different velocity dis- 
tributions would necessarily produce changes in the character of the rates of 
distortion, and (c) the external body force in the Gunaji viscometer is a con- 
stant, whereas in the Hagen-Poiseuille viscometer it is variable. 

Within the pressure and viscosity limits of these tests, the problem under 
discussion is probably more academic than practical, because the deviations 


10. Gunaji and Villemonte, “Effect of Pressure and Temperature on the 
Viscosity of Petroleum Oils”, 
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mvolved could be neglected for all practical purposes. It may well be, how- 
ver, that for higher pressures and greater viscosities, the falling-sphere 
jiscometer does not give realistic results for pipe flow problems. The im- 
jlications of this discussion are quite obvious and in line with known facts. 
for example, it might be concluded that the viscosities to be used for a given 
Ow problem should be measured by an apparatus and technique that produce 
joundary, velocity distribution, and character of distortion rate conditions 
mnalogous to those found in the problem. This conclusion is supported by the 
facts that (a) results from falling-sphere, Ostwald, rotating disk and cylinder 
vpes of viscometers cannot be used for flow in micropassages where the ad- 
sive forces become disproportionately large and entirely change the char- 
ycter oi the resistance to motion; and (b) the viscosities required for flow in 
eeve bearings seem to be some function of the rate of distortion (some- 
es referred to as “shear rate” or rate of unit strain) and can be measured 
y with viscometers which duplicate boundary, velocity distribution, and 
[stortion rate conditions. 
| The supposition that the various types of viscometry should yield different 


tesults is logical, so long as viscosity is defined as the ratio of shear stress 


b rate of distortion or strain. For a given shear stress and boundary the 
tate of distortion is a function of the molecular properties of the fluid, the 
dasic hypotheses of which have been developed by Eyring!1 and others. Vis- 
josity, therefore, is a point function and should not necessarily be a constant 
nroughout the flow regimen. The average or apparent viscosity which a par- 
cular viscometer yields can be used, therefore, only for flow situations hav- 
ng analogous flow patterns. The Stoke’s Law viscosities observed by Gunaji 
it low pressure checked the corresponding Hagen-Poiseuille Law viscosities 
iomputed from the laminar flow test section data because the Gunaji viscom- 
er was calibrated at low pressures by comparison with Ostwald viscometers 
those principle of operation is based on the Hagen-Poiseuille Law. It seems 
pgical to conclude, therefore, that when these two different types of viscom- 
ters are operated under pressures different from those at which the calibra- 
ions were made, their performances could very well be following paths which 
tre slightly divergent, for the reasons set forth above. 


7 CONCLUSIONS 

1. There was a small consistent increase in the energy loss due to pipe 
riction at 2000 psi over that at atmospheric pressure. In the absence of any 
fignificant experimental errors, it is concluded that the performance 

i. Gemant, A., “Frictional Phenomena”, Chemical Publishing Co., Inc., 
Brooklyn, N.Y., 1950. 
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characteristics of the Stoke’s Law and Hagen-Poiseuille Law viscometers 
not identical because of the different conditions in boundary, velocity distr 
tion, character of distortion rate, and external body and surface forces. 


2. It is concluded that the viscosities to be used for a given flow probl 
should be measured by a viscometer that has boundary, velocity distributi 
character of distortion rate, and external force conditions analogous to tho 
found in the problem. 


3. The standard pipe friction theory is valid for high pressures and ten 
peratures, provided the appropriate effects of pressure and temperature or 
the properties of the fluid are known. For all practical purposes, the Guna 
viscosity data (based on Stoke’s Law) can be used for pipe flow problems, 
provided the pressure and temperature ranges are reasonably close to thos 
under which the tests were made. The Appendix presents a summary of th 
data. 


The head loss in oil-hydraulics circuits for laminar flow can be estimat 
by combining the fluid property correlation equations in the Appendix with | 
Hagen-Poiseuille formula, equation 3, and for turbulent flow the Blasius fo 
mula, equation 7. The Blasius formula can be used, because the Reynolds 
Number will ordinarily be less than 150,000 for pipes that can be considers 
as “smooth”. 


APPENDIX 


A. Pressure - Density Relations for Petroleum Oils. 


The pe Bere Sous ey relations for petroleum oils are accurately repre 
sented by the Dow-Fink 2 formula 


© = Qo (I + ap - bp?) ( 


where p is the density at a gage pressure, p, in psi; Pg is the density at at: 
mospheric pressure (p=o) and any temperature; and a and b are empirical 
constants which apply for a pressure range of zero to 50,000 psi and a tem 
perature range of 20° to 220° F. The relationship betwees a and b and the 
temperature are given in Figure 10. 


B. Viscosity-Temperature-Pressure Relations for Petroleum Oils. 


The viscosity-temperature-pressure relations for petroleum oils as gin 
by Gunaji and Villemonte!3 are represented by the formula 


os e leF ees 


12. Dow and Fink, “Computation of Some Physical Properties of Lubricati 
Oils at High Pressure”. 


13. Gunaji and Villemonte, “Effect of Pressure and Temperature on the V 
cosity of Petroleum Oils”. 
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ere pl is the viscosity at the desired pressure and temperature; Uo is the 
ospheric viscosity at the desired temperature; p is the gage pressure in 
i; t is the temperature in degrees Fahrenheit; and a@,., and ¢ are correla- 
n constants which apply to a specific atmospheric viscosity at 100° F., wg. 
e relationship between the correlation constants and the specific atmos- 
ric viscosity at 100° F. is given in Figure 11. These data are applicable 
pressures up to 3000 psi and temperatures up to 180° F. 
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) SOME PHYSICAL PROBLEMS RELATED TO 
) FLOOD INSURANCE! 
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(Proc. Paper 1164) 


INTRODUCTION 
: 
he adoption of any nation-wide program of flood insurance would intensify 
need for solution of a large number of scientific and technical problems 
jompassing the fields of meteorology, hydrology, and oceanography. For- 
ately within the past two decades a large body of theoretical and experi- 
ntal knowledge has been accumulated that may be applied to obtain reason- 
> quantitative answers to some of these problems. Less fortunately, how- 
ir, all three fields are handicapped by a lack of long-time observational 
1 ae at best less than fifty years. No attempt is made in this paper 
resent more than brief descriptions of some of the more important phe- 
ena that will require the gathering of additional observations, their 
ysis and application. 
b has been assumed that the term “flood” would imply the abnormal rising 
yater of water level causing inundation of land whether by a flowing stream 
)y an ocean, sea or lake. 


Meteorology 


od-Producing Storms 


rhe severe floods in 1955 in New England and in the Pacific Coast States 
vell as those since 1950 in the Mid-western and Gulf States have empha- 
id the need for additional research in predicting rates of storm precipita- 
} produced by combinations of meteorological situations and for further 


: Discussion open until July 1, 1957. Paper 1164 is part of the copyrighted Jour- 
al of the Hydraulics Division of the American Society of Civil Engineers, Vol. 83, 
fo. HY 1, February, 1957. 


resented at a meeting of the Hydraulics Division, ASCE, October 17, 
956, Pittsburgh, Pa. 

ons. Engr., Palo Alto, Calif.; formerly associated with Parsons, 
rinckerhoff, Hall & Macdonald, Cons. Engrs., New York, N.Y. 
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statistical study of storm location with regard to ascertaining critical ge 
graphical flood-producing orientations of these meteorological variables. 

Most regions of the continental United States are subject to precipitati 
produced by storms of the cyclonic type. These storms are the result of 
mospheric waves developed along frontal surfaces separating cold and he 
polar air from warmer and lighter air originating in tropical regions. C 
clonic storms are characterized by the production of rainfall over extens 
areas; a condition met chiefly during the winter months. Upper limits of : 
clonic type rainfall rates should, as a result of extensive theoretical and 
statistical studies by the Weather Bureau, other agencies and individual it 
vestigators, be moderately reliable for application to regions east of the 
Rocky Mountains. Continued investigation will doubtless increase the reli 
bility of estimates for extreme rainfall rates for this storm type for vari 
physical and thermodynamic combinations of polar and tropical air chara 
teristics. 

Thunderstorms or convective storms are more or less common throug 
out the United States though they may locally be known by other designatic 
as for example in the arid Southwest as the “Sonora type”. They are gene 
ally responsible for the most severe floods that can occur on relatively s! 
drainage areas (500 square miles of less). The attempt to predict precis: 
where they can occur and their intensity has not been too successful. Upy 
limits of rates of precipitation and a real extent have been investigated bc 
theoretically and statistically for many regions. However, because of the 
complex internal structure much additional study will still be required be 
completely reliable estimates of their maximum flood-producing potentia. 
will be available for all areas of possible interest to a program of flood i 
surance. 

Rainfall due to the lifting action of mountain ranges on the flow of wart 
moist air over topographic barriers (orographic action) is particularly in 
portant along the West Coast. In contrast to the relatively unpredictable | 
currence of the location of cyclonic and convectional storm types, the oro 
graphic storm pattern is geographically well fixed. Also upper limits of | 
graphic storm rainfall in the Pacific Coast region can, because of rather 
initely determinable upper limits of atmospheric moisture content, be est 
mated with considerable accuracy. Reliable determinations of orographic 
storm precipitation rates for other mountainous regions of the United Sta 
are more difficult to derive. 

Various combinations of the three foregoing more of less basic storm 
rainfall types can and do occur throughout the United States. A large amo 
of investigation must still be accomplished in this difficult phase of mete: 
ological analysis to satisfy the requirements of a nation-wide flood insur: 
program. 

Large areas of the United States are subject to hurricanes, including t 
states adjacent to the Gulf of Mexico and most of the Atlantic Coastal reg 
Recent studies of rates of precipitation produced in hurricanes indicate tl 
maximum falls are at least comparable to those rates observed in severe 
thunderstorms. Possibly the rainfall of hurricanes is more intense than 
observations indicate due to the fact that high winds prevent gages from 1 
suring true rates. The relative importance of the several factors produc 
rainfall in hurricanes are not as yet well understood. Observations indic 
that the heaviest rainfall in a hurricane is produced in the right front qua 
of the storm and as a consequence for any particular watershed the flood 


; 
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oduced by a hurricane is dependent on the track of the storm. It is obvious 
at, because of the many largely undetermined inter-related meteorological 
d physical factors affecting hurricane floods, much additional data, study 
d research will yet be required to meet the technical demands of a flood 
urance program. 
The question concerning an increase of flood-producing storms in the 
ited States in recent years is a subject of much controversy among meteor- 
gists and hydrologists. At the present time no satisfactory evidence exists 
c whether a temporary or permanent change has occurred. 
| Hydrology 
The rapid development of the science of hydrology during the past thirty 
s has prepared an adequate background for the solution of most of the hy- 
logic problems that will arise in a program of flood insurance. The most 
essing need will be for more extensive observations of precipitation, rates 
snow-melt, temperature, soil moisture, streamflow, reservoir sedimenta- 
m, etc. For example, although the U. S. Geological Survey operates about 
100 gaging stations in the principal river basins there are many small 
jeams throughout the United States where potential flood damage exists but 
ere the number of flow determinations are inadequate for purposes of 
iving the technical problems arising in a program of flood insurance. A 
ilar situation exists with regard to the various types of precipitation 
asurement stations. 


| Oceanography 

| 

The great loss of life and property sustained by the coastal areas of New 

ten in September 1938 and in other locations on the Atlantic Coast since 
time has emphasized the importance of storm tides, currents and waves 
ents of destruction. High ocean water levels acting in conjunction with 
rates of flow from coastal streams can multiply property damages. along 

_ lower reaches of rivers adjacent to the ocean. The science of Oceanog- 

jhy, a relatively new field, is able to lend assistance in the solution of 

ine of these problems. 


rm Tides 

Sustained winds in the right-hand quadrants of a hurricane acting in con- 
ction with reduced atmospheric pressure near and at the “eye” of the 
rm create a rise in the level of the sea along a coast. Typical observed 


treme historical storm tides are listed below. 
) Water Level Above 


Location Date Astronomical Tide (in feet) 
v York, N.Y. September 3, 1821 13 
nnah-Charleston,S.C. August 27, 1893 Coastal islands submerged 
e Erie at Buffalo, N.Y. November 1900 8.5 
pus Christi, Texas September 14, 1919 16 
rida Keys September 3, 1935 15 to 20 
tha’s Vineyard, Mass. September 21, 1938 15 to 20 


vidence, R. I. August 31, 1954 15 
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A steady and sustained wind blowing over the surface of a relatively shal- 
low body of water creates a horizontal force on the surface of the water caus 
ing both the formation of oscillatory waves and a piling up of water. Con- 
tinental shelves along the coasts of the United States are particularly sus- 
ceptible to this piling up of water under the action of onshore winds. In addi- 
tion variations in atmospheric pressure cause water levels to rise when the 
barometric pressure is low. This latter action is magnified if the barometri 
low travels at a speed approaching that of the propogation speed of the water 
wave (resonance effect). The combined result of the wind tide and the baro- 
metric tide is the total storm tide. 

Theoretical procedures are available for computing storm tides that 
should permit reasonable approximations of maximum possible water levels 
for any coastal location. For estuaries the total rise at any given point would 
be a combination of static tide and the tide due to seiches (oscillation of the 
body of water in the estuary). 

Studies indicate that, with the possible exception of the coast of Southern 
New England, no portions of the coastal waters along the Atlantic seaboard 
and Gulf of Mexico have as yet experienced the greatest possible storm tides 
due to hurricanes. 


Currents 


Under normal circumstances reversing currents exist at the mouths of all 
estuaries. These currents extend their influence along the coasts on both 
sides of the opening to the sea. 

During severe storms creating offshore or onshore winds there is a tem- 
porary variation of water surface due to the transport of surface water in the 
ocean and in an estuary. Frequently during violent storms the additional 
volume of water transported in and out of an estuary generates exceptionally 
strong reversing currents that are hazardous to navigation and can occasion- 
ally cause marked shoreline changes due to erosion and transport of beach 
material. 

For many locations data concerning the magnitude of extreme current 
velocities created by storms is entirely lacking, a deficiency that might re- 
quire rectification depending on whether the resulting damage is classified 
as due to floods. 


Waves 


Storm-generated waves superimposed on high static water levels caused 
by storm tides can be agents producing very severe damage to coastal im- 
provements of all types. Forces exerted by breaking waves in particular are 
extremely large and until fairly recently the great magnitude of the pressures 
exerted on structures has often not been adequately appreciated. Past and 
continuing observation and research on the complex properties and action of 
waves on structures gives the engineer adequate tools for safe design of 
coastal structures. 


Tsunami 


Although the majority of water waves are the result of meteorological | 
causes, mention should be made of the wave produced by submarine earth- 
quakes and designated “tsunami”. In the open sea these seismic sea waves 
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have small amplitude and travel at high speed. On approaching mildly sloping 
shores the waves rapidly increase in amplitude and are very destructive. In 
bays, harbors and estuaries the “tsunami” may initiate seiches whose periods 
of oscillation depend on the depth and shape of the body of water. 

Whether the seismic sea wave is a serious threat to coastal areas of the 
United States is not definitely known from observation, due to the shortness of 
tidal records. It may be a subject worthy of investigation if “tsunami” should 
be classified as a flood phenomenon in any insurance program. 


SUMMARY 


EE 


An attempt has been made to outline in general terms some of the techni- 
cal and scientific problems that would be met in a nationwide flood insurance 
program. Additional observational data related to meteorology, hydrology 
and oceanography would be a basic requirement. Although a large body of 
: knowledge now exists that is readily available to assist in the immediate 
: solution of some of the complex problems, an expanded program of investiga- 
tion and research would be a necessity. 
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TECHNICAL PROBLEMS OF FLOOD INSURANCE* 


H. Alden Foster! 
(Proc. Paper 1165) 


ABSTRACT 


This paper discusses Flood Insurance from an engineering standpoint with 
particular emphasis on determination of flood probability from relatively 
short records, estimates of mean annual flood damage to particular proper- 
ties for establishing self-supporting premium rates, spreading the risk of 
flood damage over a large number of policies, and related technical problems 
involved in setting up any program of flood insurance. 


INTRODUCTION 


The problem of bringing financial relief to property owners from the ef- 
fects of damage by floods has been receiving increasing attention in recent 
years. The private insurance companies as well as Congress have given 
much study to this question. In much of the discussion that has taken place 
on this subject, it has been apparent that there has been considerable misun- 
derstanding of the technical problems involved and the difficulties that would 
have to be overcome in establishing any program of flood insurance or in- 
demnity, either on a private or government-sponsored basis. 

The speaker proposes, therefore, to discuss this subject in a general way 
and to bring out the basic problems involved, without attempting to prove 
whether flood insurance is actually feasible. This latter question will be dis- 
cussed by another speaker. ; 


Note: Discussion open until July 1,1957. Paper 1165 is part of the eonynlanes Jone 
of the Hydraulics Division of the American Society of Civil Engineers, Vol. 83, No. 


HY 1, February, 1957. 

* Presented at a meeting of the Hydraulics Division, ASCE, October 17, 1956, 
Pittsburgh, Pa. This paper is based largely on Proceedings Separate No. 
483, “Flood Insurance,” by H. Alden Foster, August, 1954. 

1. Princ. Associate, Parsons, Brinckerhoff, Hall & MacDonald, New York, 
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Basic Problem of Flood Insurance 


When an insurance company issues a policy to idemnify the policy-holder 
against a certain type of loss, such as by fire, the policy provides that the 
owner will pay a definite annual premium to the company as compensation for 
the risk carried by the latter. This premium involves two parts: (1) the 
basic or “pure premium” which represents the average annual loss sustained 
by the company, and (2) a “loading charge” intended to cover the company’s 
cost of doing business. 

The basic problem in flood insurance, therefore, is to determine the Aver- 
age Annual Damage that may be expected to occur to each property on which 
insurance is issued. This problem may be considered in three steps: 


1. Estimating the average annual frequency or probability of occurrence 
of floods of a given magnitude at the locality under consideration. 


2. Estimating the dollar value of damages caused to each property by a 
flood of any given intensity. 


3. Determination of the Average Annual Flood Damage to the property, by 
combining the results of the above estimates. 


Flood Frequency 


Statistical Methods 


Estimates of future flood occurrences must be based on past experience, 
modified by the anticipated effects of any changes in controlling conditions © 
that may be expected to occur in the future. While it is impossible to fore- 
cast the actual flood intensities that will be experienced in any future year, it 
has been possible to make estimates of the probability of occurrence of floods 
of various intensities, by statistical analysis of past records. 

Statistical methods of considerable refinement have been used for many 
years in various fields, particularly in life insurance where technical pro- 
cedures of considerable value have been developed. Somewhat similar 
methods were introduced into the field of hydrology about forty years ago, 
and have received intensive study by many engineers within the past twenty 
years. However, there is a fundamental difference between the statistical 
methods used by engineers and those that have received general acceptance 
by statisticians. 

In most statistical work, the data available for analysis run into a very 
large number of items or individual events, as for instance with life insur- 
ance. But in hydrology and particularly in the study of flood records where 
the maximum flood occurring each year is considered, the number of items 
available in any record is quite limited, seldom extending over more than 
fifty years and often being limited to twenty years or less. 

Statistical methods, at best, must necessarily yield approximate results. 
Where the available data are limited in quantity, the degree of precision in 
the calculated results will become progressively smaller. This is an impor- 
tant consideration in evaluating any program of flood insurance. 


Selecting the Curve-of-Best-Fit 
Various methods have been advocated for plotting the flood-probability 
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curve. The purpose of all of these methods is to establish an approximate 
relationship between the magnitude of any particular flood and the probability 
of its occurrence in any future year. Some of the methods that have been ad- 
vocated for this purpose are adaptations of methods extensively used in 
general statistical analysis; others are based on formulas developed specifi- 
cally for use in hydrological investigations. It is beyond the scope of the 
present discussion to make a detailed comparison of these methods. 

When the flood data are plotted on a diagram showing flood magnitude 
versus probability of occurrence, it is desirable to replace the individual 
items by a smooth line called the “curve-of-best-fit.” While the advocates 
of some of the plotting methods maintain that their curves have a reliable 
mathematical basis for determination of this curve, it appears doubtful 
whether these claims can be fully substantiated. In the opinion of the 
speaker, the best proof of a scientific basis for a given probability curve is 
that it gives a reasonably good fit to the recorded data. The practical tests 
of reliability are: (1) which method gives the best representation of actual 
data; and (2) which method is most convenient to use. 

The various proposed methods of curve-fitting generally produce reason- 
ably good results with a given flood record, when used to express probabili- 
ties or recurrence-intervals within the length of the record. For example, 

a 20-year record could be used to obtain a reliable estimate of 10-year (10% 
probable) flood, regardless of which method is used for obtaining the proba- 
bility curve. But when an attempt is made to extrapolate the flood-probability 
curve much beyond the recurrence interval covered by the record, consider- 
able differences are revealed in the several methods. However, these differ- 
ences are far outweighed in magnitude by the so-called “errors of sampling.” 


Errors of Sampling 


In the language of statistics, a record of annual floods at a particular loca- 
tion is a “sample” of all floods occurring at that location over a long period 
of years, the long-time record being known as the “population” and the indi- 
vidual recorded floods being items of the population. The items of the re- 
corded sample may be analyzed by the various probability methods previous- 
ly mentioned, and a smooth probability curve may be established as giving 


_ the best representation of the sample. But there is no assurance that a sam- 


ple record obtained in the future and covering the same length of time will 
produce the same probability curve. In a program of flood insurance, it is 
important to determine the extent to which the future record may differ from 
the available sample. 

Methods have been developed in the Theory of Statistics whereby the prob- 
able errors of future samples, or the “spread” of the items in the future 
sample as compared with the available record, may be reasonably estimated: 
such methods show that the relative magnitude of the errors decreases rapid- 
ly as the number of items in the record increases. For records of less than 
twenty years, the possible errors may be sufficiently large to appreciably af- 
fect any estimates of mean annual flood damages. : : 

A conspicuous example of this situation is found in the “Hurricane Floods 
of August 1955 in the northeastern states. These floods in many locations 


greatly exceeded any that had occurred within the previous century. A prob- 


ability curve obtained from a flood record in such locations prior to 1955 : 
would be substantially different in shape from one based on a record including 


the 1955 floods. 
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Changes in River Conditions 


Where regulating works have been constructed on any river and have the 
effect of reducing the peak discharges of floods, the probability of occurrence 
of floods on the river will be affected. Any studies of flood-probability on 
such a river based on records prior to the construction of the regulating 
works would have to be revised before they would be useful in a subsequent 
program of flood-insurance. Obviously, the amount of detailed study re- 
quired for preparation of flood-probability estimates at a given location will 
be appreciably increased if the discharge of the river is affected by flood- 
control structures further upstream. Similarly, changes in watershed condi- 
tions, either in the past or anticipated in the future, which would cause in- 
creased flood discharge, such as removal of forests or the conversion of land 
to urban uses, must be given consideration in the determination of flood 
probabilities. 


Conclusions as to Methods of Curve- Fitting 


After studying and comparing various methods of curve-fitting that have 
been proposed for hydrological studies, the speaker has not been able to 
reach any definite conclusions regarding which one is the best for use in 
flood-insurance studies. Within the time limits of the actual record, the 
several methods seem to give substantially equivalent results. However, 
when the curves are extrapolated beyond the highest floods of record, there 
is considerable variation among the results obtained. 

Since one of the principal objects in preparing a flood-probability curve 
for insurance purposes is to permit the extrapolation of flood frequencies be- 
yond the time limits of the record, the question of the amount of possible er- 
ror in the “sample” for less frequent floods is of particular importance. 

Such errors may become relatively large, especially when the record covers 
only a few years. In fact, these errors of sampling far outweigh in magnitude 
any differences between the several types of theoretical probability curves 
that have been mentioned. Hence the selection of the most suitable method 
comes down to a question of convenience in application. Fortunately, a pre- 
cise determination of the magnitude of floods with very low probabilities (the 
floods with very long recurrence intervals) is not of great significance for 
computation of insurance premiums, since these rare floods have only a 
slight influence on the value of the “average annual flood.” 


Flood Damages 


The second part of the problem of estimating average annual flood losses 
for a particular property is the computation of the probable damages to the 
property in question that would be caused by a flood of a given intensity. This 
matter must be considered in two parts: 


1. To what depth will the property be inundated by the flood? 


2. How much damage will be caused to the property when the water 
reaches that depth? 


The first question involves the relation between “flood stage,” or flood 
elevation, and flood intensity; the second involves the relation between flood — 
stage and damages for the property in question. 
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Flood Stages as Related to Particular Properties : 


The work required for the preparation of diagrams of flood stage at a par- 
ticular property, as related to the intensity of the flood, involves; 


1. Preparation of necessary topographic maps showing the location and 
elevation of all properties under consideration at a relatively large 
scale, 


2. Preparation of River Stage-Discharge curves, showing the elevation of 
the water surface at the gaging station for various intensities of dis- 
charge. 


The preparation of the necessary maps and stage-discharge curves for a 
given locality may require a great amount of field surveying as well as office 
preparation, and will involve considerable time and expense for preparation 
of the insurance program. This must be incurred regardless of the amount 
of flood insurance subsequently issued in the area. 


3. Determination of the MSL elevation of each property under considera- 
tion, It is important to determine the “elevation of zero damage,” or 
the flood stage at the property that must be reached before any materi- 
al damage will occur. 


Classification of Flood Damages 
Damages to property are classified as follows: 


1. Direct losses, consisting of physical damages to property and goods, 
measured by present day cost of repair or the replacement in kind, and 
the cost of cleanup and for moving goods. The direct losses may also 
be subdivided according to the type of use for the structures,—such as 
residential, commercial, industrial, etc. 


Indirect losses, consisting of the value of services or uses lost or made 
necessary by reason of flood conditions, not chargeable to direct loss. 
These include losses of business and wages, costs of relief and similar 
losses, both within and without the flood area, during the period of the 
flood and subsequent rehabilitation. It seems doubtful whether indirect 
losses of these types would be included in a program of flood insurance. 
On the other hand, indirect flood losses have commonly been included in 
determination of the benefits produced by construction of flood-control 
projects. 

3. Depreciation losses have occasionally been included as a credit for the 
economic justification of flood-protection works. The total depreciation 
of a property that has been flooded equals the value of the property be- 
fore the flood minus the value after the flood. It appears doubtful, how- 
ever, whether such depreciation losses would need to be considered in 
connection with flood insurance if physical damage to the property has 
been compensated for under “direct damages.” 


pe 


Estimates of Damages 


The determination of the amount of damage that would be caused to the 
property by a particular flood stage is a matter of appraisal. The most ead 
able estimated of damage will be for floods of record from which the owner o 
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the property has definite records of losses sustained. Generally the apprais- 
er must make his own estimate, which may involve a detailed examination of 
the property and the installed equipment, followed by estimates of losses to 
be sustained by each piece of machinery or equipment, and of the cost of 
cleaning up. 

Damage to buildings generally will be less than the total replacement cost 
unless the type of construction is such as to result in total destruction, such 
as by undermining the foundation. 


Dynamic Effects of Floods 


The stage-damage estimates previously described apply particularly to the 
inundation type of flood, where the water surface gradually rises until the 
property is partially or fully submerged. This type of flooding will occur 
where the slope of the flood-profile is moderate. However, where the surface 
profile of the river is relatively steep, or where the flood waters are tem- 
porarily held back by obstructions, high velocity currents may be induced 
that would cause serious damages. A simple inundation may cause only par- 
tial damage to a structure, whereas the dynamic effect of rapidly flowing 
water may bring about its complete destruction. These effects are so diffi- 
cult to estimate in advance that they can only be approximated by reference 
to past experiences in a given locality. 

The type of flood damage may be greatly influenced on a particular stream 
by the character of the storm producing the flood, as was conspicuously 
demonstrated by the “hurricane floods” of August 1955 in the northeastern 
states. On the larger rivers, such as the Delaware, the damage was caused 
largely by water transmitted from upstream at greater rates than could be 
handled by the natural river channel, resulting in overflow of the river banks 
and inundation of the adjacent flood plains. On the smaller streams there 
was apparently extensive sheet runoff caused by unprecedentedly high rates 
of rainfall resulting in “overland floods,” with widespread damage and loss of 
life in sections that had never previously been subject to severe floods within 
the period of record. 


Average Annual Flood Damage 


The average annual flood damage experienced by the owner of a particular 
property is the amount which the insurer would have to collect each year 
merely to pay the losses that will occur during a period of years. Assuming 
that the property is restored after each flood experience, the total losses 
over a period of years divided by the number of years would be considered 
the average annual damage. The calculation of this quantity is made as fol- 
lows: 

The diagram of flood stage vs. flood discharge at the gaging station cor- 
responding to the given property (also called the “rating curve”), is combined 
with the diagram of flood discharge vs. probability, for the same gaging sta- 
tion, to obtain a diagram of floodstage vs. probability (Fig. 1). This is ap- 
plicable to all properties in the local area. 

The latter curve is then combined with the diagram of floodstage vs. 
property damage previously described to obtain a curve of flood-damage vs. 
probability for the particular property (Fig. 2). This curve shows the proba- 
bility of occurrence of flood damages equal to, or exceeding, various 
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amounts. It can be readily shown that, if the curve is plotted on a numerical 
scale, the total area under the damage-probability curve is a measure of the 
average annual damage caused by all floods which may be anticipated over a 
long period of years. 

The value of the Average Annual Flood Damage so determined must neces- 
sarily reflect the approximations involved in preparing the several diagrams 
of flood stage, probability and damage. The least reliable of these is likely 
to be the diagram of stage vs. probability. 


Test Calculations of Average Annual Flood Damage 


The procedures previously described have been used to calculate the aver- 
age annual flood damages to a considerable number of properties at Lowell, 
Mass. This location was selected because much of the basic data were avail- 
able in the files of the U. S. Corps of Engineers. The ratio between the mean 
annual loss for each property and the maximum probable loss for the proper- 
ty was used as a basis of comparison, since it would be a factor in establish- 
ing the premium rate for a flood-insurance policy. Out of 54 industrial prop- 
erties so investigated, the value of this ratio varied from 0.4% to a maximum 
of 17.5% with an average value of 4.2%. Similar results have been obtained 
from studies in other localities. This indicates that the cost of flood insur- 
ance to a property owner would be high if such a program were feasible. 

Another important consideration is the great variation in the value of the 
ratio of mean annual to maximum flood losses, which would either require a 
great variation in the premium rates charged on different properties in the 
Same area or would force the more favorably located properties to carry a 
considerable share of the insurance cost on the properties subject to greater 
risk. 


Financial Reserves 


Basing the insurance premium on average annual flood damage makes it 
theoretically possible for the insurance company to reimburse itself for 
losses over a long period of years. But this does not eliminate the possibility 
that the maximum possible flood may occur in any year; and unless the under- 
writer has set up reserve funds sufficient to pay such maximum loss at any 
time on all insured properties subject to flooding by a particular river, the 
financial stability of the insurance program would be jeopardized. Conse- 
quently, the total reserves required to assure the solvency of the underwriter 
for any general area would be equal to the sum of the maximum insured 
losses of all the insured property in that area. 


Spreading the Risk 


Obviously, there is no way by which the total risk may be reduced by 
spreading it over a large number of properties when they are all located in 
the same general flood area, since a major flood will cause damage to all in- 
sured properties located within the area affected by the river in question. In 
this respect, flood insurance differs materially from fire insurance. Under 
normal conditions, only a small proportion of insured properties in a particu- 
lar city suffer damage by fire in any one year, SO that the risk of loss can be 
spread over a large number of policies. It is only in the case of a 
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catastrophe, such as the Chicago fire of 1871, that a large number of proper- 
ties are damaged by fire at one time, whereas every major flood may be con- 
sidered a catastrophe since it would result in damage claims from all insured 
owners within the affected area. 

If the insurance program were placed in operation in widely separated 
sections of the country, it is probable that some beneficial effect of “spread- 
ing the risk” may be realized. This is because major floods do not occur in 
all parts of the country in the same year. A very approximate study of flood 
loss statistics for the entire United States indicates that the reserve set 
aside to cover major flood losses in any one drainage basin might be reduced 
by two-thirds if the insurance program covered the entire country and was 
not limited to individual basins. 

It is important to note that this question of spreading the risk applies only 
to the maximum flood losses and does not affect the mean annual loss for any 
property. The total flood loss for the entire country is not reduced in any 
year, even if the insurance program covers the whole country. The practical 
effect of a widespread coverage by insurance, therefore, would be only to re- 
duce the amount of reserves required to cover large losses in particular 
years, 


Effect of Flood Forecasting 


Forecasting of floods on the more important rivers has been carried out 
by the U. S. Weather Bureau for many years, the extent and reliability of the 
forecasts having appreciably increased in recent years. If flood insurance 
were to be established in this country, it would doubtless be greatly benefited 
by an efficient forecasting service, since damages could be reduced in many 
cases by emergency measures taken by the property owners and based upon 
advanced warnings supplied by the flood-forecasting service. However, the 
following questions would remain to be answered: 


1. If the property owner fails to take due advantage of the forecast by de- 
laying to take the necessary precautions or neglecting them entirely, 
can the insurance company refuse to indemnify him for part or all of 
his losses? 

2. How can the insurance company be sure that the flood forecast will be 
reliable or will be issued in time so that precautionary measures can 
be taken by the covered property owner? 


Special Problems 


Numerous special problems would have to be solved in setting up a pro- 
gram of flood insurance, a few of which may be mentioned: 


Types of Floods 


Floods are produced by numerous causes and may be classified into vari- 
ous types. It may be difficult in some cases to determine where damage to 
property has been caused by direct rainfall or by some type of flood. For _ 
example, certain areas that normally would not be considered as located ina 
flood-plain may be subject to damage from heavy rainfall flowing over steep 
ground slopes in the form of “sheet runoff.” There might be some doubt 
whether this condition should be properly classified as a “flood,” or whether 
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tis a case of direct damage by rainfall. Similarly, the “mud flows” or dis- 
harges of large quantities of debris carried out on the alluvial fans at the 

ouths of steep canyons, such as in the Los Angeles area of California, may 
Pause widespread damage to structures, even though the discharge of water 
[a not be sufficient to cause great damage by itself. 

Some types of flood are difficult or impossible to forecast from a proba- 
vility standpoint. In a flood insurance program, it would doubtless be desir- 
able to place certain limitations on the type of flood covered by the policy 
\Ssued on any property, although practical considerations might make it diffi- 
sult or impossible to impose such limitations. 


Effect of Insurance on Land Use 


The owner of a property covered by flood insurance may be inclined to in- 
srease its use on the basis that, if a flood comes, the additional loss to which 
1€ would be subjected would be compensated to him. This condition might be 
aken care of by making an annual inspection of the property and, if neces- 
sary, increasing the premium charge to correspond with the increased risk. 

Uneconomical development of land subject to periodic inundation may be 
-ontrolled to some extent by the fact that insurance premiums on such prop- 
srties would necessarily have to be relatively high, thereby having a tendency 
0 eliminate properties of that type from the insurance program. 


tifect of Protective Works 


Reduction of the risk of damage by flood to properties on rivers where 
lood-protection works have been constructed will depend to a considerable 
‘xtent on how the works are operated during a particular flood. While the 
yvorks may have been designed on certain assumptions as to the manner of 
heir operation, there is no guarantee that they always will be operated in 
onformance therewith. If insurance is issued on the basis of those assump- 
ions, the additional risk involved in possible changes in the operating pro- 
Tam would have to be assumed. Similar additional risk may have to be car- 
ied in the case of lands protected by levees that are not properly maintained. 

Owners of property located downstream from flood control works may feel 
false sense of security by assuming that the degree of protection from 
hich they benefit is greater’ than that actually provided in the design of the 
forks. This may result in an increased use of land that is still unprotected, 
nd greater risk of damage from floods. 


Other Applications of Probability Methods 


The general methods outlined herein for determining the Average Annual 
'lood Damage have been extensively used to determine the Average Annual 
nefits accruing to a flood control project as a result of protection furnished 
) various properties in the flood basin. The purpose of such estimates is to 
istify construction of the project on an economic basis, estimating future 
enefits by reference to past flood experiences. If the indicated benefits ex- 
eed the annual cost of the project, the latter is considered to be justified in 
ecordance with law. If, in the future, greater floods occur than had been 
nticipated from past records, the project may still have been justified since, 
ithout the protective works, the average annual damages would doubtless 
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have been considerably greater. The government agency that constructs the 
protection works is not held liable for losses sustained by property owners 

if a flood occurs exceeding in magnitude that for which the project was con- 
structed. 

With flood insurance, on the other hand, the insuring agency agrees to re- 
imburse the property owner for all losses that may occur in the future, unde 
the terms of the policy. If the resulting average annual loss exceeds that 
estimated from past flood records, the insurer necessarily takes a loss. 
Consequently, the question of the reliability of probability estimates of futur: 
floods is of much greater significance in a program of flood insurance than i 
justifying a flood protection project. 


SUMMARY 


If a program of flood insurance were set up, the following factors require 
careful consideration, due to their influence on the charges required for pro 
tection on individual properties: 


1. The Average Annual Loss or Damage, which is the amount of money 
that would have to be set aside from premium charges each year to 
cover losses chargeable to the policy. It equals the sum of all losses 
expected to occur during a considerable number of years in the future 
divided by the number of years, and may be estimated by a study of the 
probability of occurrence of floods of various magnitudes combined wi 
estimates of the amount of damage which the various floods would cau; 
to the property. 


2. The Errors of Sampling, which are a measure of the extent to which 
future flood records may vary from the “sample record” on which the 
estimate of Average Annual Flood Damage is based. No precise esti- 
mate of future floods can be made in advance; and the shorter the avai 
able record the greater the possibility of variation in the future. 


oo 
. 


Sufficient Financial Reserves must be available to the insuring agency 
to equal the sum of all the maximum damage items on insured proper- 
ties located in the same general flood area. If the insurance program 
covers a sufficiently large portion of the country so that maximum 
flood losses will not occur in all portions of the insured areas in the 
same year, the amount of reserve funds can be reduced to some exten 
Such reduction, however, would not involve any change in the average 
annual flood losses for which provision must be made in establishing 
the premium rates. 


There is no constant ratio between average annual flood loss and the 
amount of insurance on a property, even for similar types of properties. 
This ratio is an important factor in establishing the annual charge, and will 
vary greatly with the location of the property with respect to the river chan- 
nel. In general, premium charges will have to be determined individually fc 
the various properties. Rates for flood insurance will have to be much > 
greater than for other types of insurance. | 

There are numerous other problems involved in setting up a program of 
flood insurance, the effects of which are difficult to assess in advance be- | 
cause of lack of actual experience in this field of insurance. 
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SYNOPSIS 


The basic principles of frequency analysis of streamflow data are re- 
yiewed pointing toward the advantages and applications of these simpler 
orms of statistical analysis. Computational procedures are only briefly out- 
ined, but a list of references indicates where details can be found. Fre- 
juency analysis of high flows (floods), low flows (droughts), and all flows are 
lescribed. A group of figures illustrates methods of expressing and inter- 
reting streamflow data subjected to frequency analysis. 


INTRODUCTION 


Few civil engineers profess to be statisticans; in fact, few claim even an 
lementary knowledge of statistics. Lacking formal instruction in a subject, 
t is natural to avoid the unknown. Yet, an awareness of the advantages of 
tatistical analysis encourages the application of such methods. Too often 
nnecessarily complex terminology and a variety of means of obtaining simi- 
ar results prove to be the most confusing elements. Thus, a presentation of 
asic methods, omitting specific details and derivations but emphasizing ap- 
lications, may be worthwhile. 

An example of the application of simpler forms of statistics is that re- 
uired for frequency analysis of streamflow data. If data are available, the 
ivil engineer is faced with analyzing a large collection of numbers—usually 
1ean daily discharges—to obtain information related to an engineering prob- 
2m. It may involve a structure to store, pass, or divert water, or it may 
oncern water supplies for consumption or power generation. In any event, 
1e problem reduces to determining the frequency distribution of streamflow 
ata. Frequency may be defined as the number of flows of a particular mag- 
itude as compared to the total number of all such flows, whereas a frequency 
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distribution is an arrangement by magnitude of frequencies. For example, in 
a given 100-year record there are 100 annual peak flows, or floods. The 
smallest of these flows is equalled or exceeded every year and therefore can 
be designated the 1-year flood. Similarly, the fifth largest flood occurs on 
the average every 20 years and becomes the 20-year flood; the highest flow, 
the 100-year flood. These recurrence intervals Tr may be defined by the 
familiar equation 


jt (1) 
m 

where m is the order number when the values are arranged in decreasing or- 

der and n is the total number of events. It should be noted that no chronologi- 

cal sequence is implied but only an average interval of occurrence. 

Frequency distributions can be converted to frequency curves by simple 
arithmetic procedures. Many plotting schemes have been devised to facilitate 
interpretation of frequency curves; familiarity with these is not as important 
as understanding the limitations of any one method. 

The following three sections describe the applicability of frequency analy- 
sis to high flows, or floods; to low flows, or droughts; and to all streamflows. 
References at the end of the paper provide an introduction to the extensive 
literature available on the subject. 


Frequency Analysis of Flood Flows 


Flood frequencies are developed for design of culverts, bridge waterways, 
and spillways, for preliminary planning, and for economic studies. Available 
data usually consist of streamflow records published by the U. S. Geological 
Survey in annual Water-Supply Papers as mean daily discharges and annual 
peak flows. Floods may be defined either as peak or volume floods depending 
‘upon the purpose of the analysis. For floods occurring above dams having a 
relatively large storage capacity, the frequency of peak flows is incidental to 
that of flood volumes; for small dams containing minimal storage, the con- 
verse is true. A peak flood is simply the highest instantaneous discharge; 
whereas, a volume flood is that high rate of discharge persisting for a speci- 
fied duration, such as a 1-day, 3-day or 10-day flood. Volume floods are ob- 
tained by accumulating mean daily discharges over the time interval and con- 
verting to acre-feet. - 

To develop a peak flood frequency curve, peak floods are selected from a 
streamflow record. Ordinarily a minimum of 10 to 20 years of record are 
necessary to approximate a frequency curve. Floods can be selected either 
by the highest flood in each water year (the annual flood method) or by all 
peak floods above the lowest annual flood (the basic flood method). I peak 
floods by both methods are arranged by order of magnitude, an array such as 
that shown in Fig. la for a 20-year record(7) will be obtained. The 20 annua’ 
floods are indicated by dots and the 20 highest floods known as the annual ex- 
ceedance floods, are designated by horizontal bars. For comparison these _ 
two sets of selected floods are replotted in Fig. 11. It can be seen that the 
annual exceedances form a population closely approximating the annual maxi 
ma for the higher floods but deviating positively as the lower floods are _ 
reached. It can be shown from probability theory(7,22) that the recurrence 
interval obtained by the two methods differs most from the lowest floods and 
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least for the highest floods. Fig. 2 shows the relationship of the two methods, 

graphically. For a 10-year flood the recurrence interval discrepancy 

amounts to 5 percent and reduces to 1 percent for a 50-year flood. 

For temporary or small structures where designs are based upon flood 

_ frequencies of about 10 years or less, the basic flood method should be 

_used;(26) for large permanent structures designed for larger floods either 

| method is satisfactory but the annual flood method is preferable because it is 

shorter and avoids the problem of establishing the independence of consecu- 
tive floods. 

| A frequency distribution for a series of floods, grouped by classes to form 

a bar graph, is shown in Fig. 3a. If the number of items were increased and 
the class size reduced, the theoretical curve would be formed. The unsym- 
metric shape is typical of flood frequency distributions. Now, if these values 
are accumulated from the highest to the lowest and divided by the total num- 
ber of flood values, a cumulative frequency curve (commonly referred to 
simply as a frequency curve) is formed (Fig. 3b). The ordinate expresses 
the probability of recurrence of a flood equal or greater than that specified 
by the abscissa. For ease in reading and extrapolating such curves, special 
coordinate scales have been developed to transform the usual S-shaped curve 
nearer to a straight line (Fig. 3c). The particular theoretical frequency 
curve employed determines the type of scale required to produce a straight 
line. If the observed flood frequency curve differs from the theoretical one, 
a curved line will result. 

A common problem upon developing a flood frequency curve is the need to 
estimate the magnitude of floods having recurrence intervals greater than 
that of the length of record. The extreme curvature displayed by most fre- 
quency curves plotted on rectangular coordinate paper makes extrapolation 
difficult; therefore plotting the curve on probability, log-probability, semi- 
log, or log-log paper is sometimes helpful in reducing the curvature. Most 
extrapolation by eye of frequencies is arbitrary, hence it does not provide a 
sound basis for estimating frequencies of high floods. To overcome this dif- 
ficulty analytical techniques have been developed for fitting a cumulative fre- 
quency curve to specified types of frequency distributions. The contributions 
in this connection of men such as Foster, Hazen, Powell, Goodrich, Slade, 
Potter, Gumbel, Chow, and Geyer are well-known to many engineers. Essen- 
tially, these efforts have produced a variety of ways for extending a flood 
frequency curve from a limited length of record. Computational procedures 
in most cases are straight-forward—tabular forms reduce the work in many 
instances to almost a mechanical procedure. 

With the aid of statistical parameters determined from the streamflow 
data and special plotting papers, it is possible to estimate say the 10,000- 
year flood from a record of 20 years duration. That such practices are 
dangerous has been well emphasized by the ASCE Subcommittee on Review of 
Flood Frequency Methods.(26) A short streamflow record obtained during a 
climatically atypical series of years is a poor sample to use for estimating 
large floods of infrequent occurrence. Generally speaking, the smaller the 
frequency required, the more accurate the estimate becomes. t 

Chow(6) compared the equations for several of the better known analytical 
procedures developed for extending frequency curves. He showed that differ - 
ent methods based on identical data will produce different flood magnitudes 
for a given frequency. Nevertheless, he went on to demonstrate that the dif- 
ferent methods are similar in structure; that is, a general formula was 
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applicable to all of the methods with a “frequency factor,” contained in the 
formula, constituting the only variation. 

A flood frequency curve for a Kansas river(19) is shown in Fig. 4. This 
curve is plotted on Gumbel probability paper(25) and was computed from the 
recurrence intervals given by Eqn. 1. If the observed annual floods fitted the 
theory of extreme values, which is the basis for the Gumbel frequency dis- 
tribution of hydrologic data, they would fall along a straight line. It can be 
seen that for practical purposes this does occur except for the highest flood 
of record which plots noticeably above the fitted straight line. Such devia- 
tions of one or more of the highest floods are quite common and, because of 
their confusing position, have called forth a variety of explanations for their 
existence. Depending upon the basin in question, one or more of three situa- 
tions may apply; (a) the hydraulics of large floods, involving channel charac- 
teristics and valley storage, may differ from that of small floods; (b) the 
cause of the high floods may be different from that of low floods; the con- 
trasting nature of rain floods and snow floods in mountainous basins furnishes 
a good example; and (c) based on the laws of probability, an unusually large 
flood can occur in a relatively short length of record. Causes (a) and (b) 
represent changes within the annual flood record at a given gaging station and 
suggest that two different populations are intermingled; it follows, therefore, 
that not one but two straight lines of differing slope should be used to fit such 
flood frequency curves. On the other hand, cause (c) implies that if the 
record were sufficiently long the unusually high flood would plot on the 
straight line with a much longer recurrence interval. Ignoring such floods is 
then an obvious recourse on the frequency curve. Yet, it is well to note that 
just as a 100-year flood will be outstanding in a 20-year record, so will a 
1000-year flood in a 100-year record. 

The recurrence interval for a given flood magnitude was defined as the 
average time interval within which the flood will be equaled or exceeded once. 
Thus, a common assumption is that in a 100-year record it is highly probable 
that one flood will occur representative of a 100-year recurrence interval. 
But, from probability theory, it can be stated that the chances of this occur- 
ring are small. To illustrate this, data from the ASCE flood frequency re- 
port(26) are plotted in Fig. 5. Here, a mythical 10,000-year record was as- 
sumed which should contain 100 floods having a recurrence interval of 100 
years. Fig. 5 shows that only 37 centuries contained one 100-year flood, 
whereas another 37 had none and 26 others had more than one. 

In hydrologically homogeneous areas frequency curves developed on gaged 
streams can be employed to estimate flood frequencies on ungaged areas. (26) 
Flood magnitudes are plotted against drainage area for equal frequencies. 
Lines fitting the data form a family of curves for different frequencies, there- 
by enabling frequency curves to be estimated for an ungaged basin within the 
area. Fig. 6 illustrates the result for the Atlanta, Georgia, metropolitan 
area.(5) Flood volumes can be treated in a similar manner; Fig. 7 shows 
the relation among discharge, duration, frequency, and area on one set of 
coordinates for a California stream. (4) 

When basins in nonhomogeneous areas are compared on a discharge-area- 
frequency plot, a scatter develops indicating that parameters other than area 
are involved. Topographic, geologic, and climatic factors cause these differ- 
ences. For example, 25-year floods in Kansas,(19) shown in Fig. 8, display 
definite regional patterns. If the particular factors and their interrelations 
were understood, it would be possible to extend the applicability of flood 
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frequency estimates on a more general basis. Research has indicated that a 
better understanding of this important problem is possible. Efforts now un- 
derway by the Soil Conservation Service, Geological Survey, and Bureau of 
Public Roads hold promise of better estimates of flood frequencies in the 
future. 

An important aid for determining the reliability of a theoretical flood fre- 
quency curve is furnished by confidence limits. Gumbel(16) developed the 
technique for finding confidence limits of annual floods after which Chow(7) 
reduced the computations to a simple procedure. Confidence limits are lines 
constructed above and below a frequency curve such that approximately 68 


percent (one standard deviation) of all observed floods will fall within the 


band. The confidence limits are functions of the flood frequency, the length 
of record, and the standard deviation of the observed flood frequency distri- 
bution. For illustration, the control curves defining the confidence limits for 
a peak flood frequency curve of a Utah stream are shown in Fig. 9. These 
indicate that from the available data there is a 68 percent probability of esti- 
mating the 100-year flood within an error of + 20 percent. 


Frequency Analysis of Low Flows 


Low flows can be analyzed on a frequency basis in an identical manner of 
that of flood flows. Instead of selecting peak flows, the lowest recorded dis- 
charges each year are employed. If droughts are defined as the annual mini- 
ma of discharges,(17) drought frequency curves are obtained. Then too, 
droughts of varying durations correspond to flood volumes. Information of 
this type has important engineering applications for water supplies dependent 
upon stream-flow, for hydroelectric power developments, and for estimating 
storage requirements. 

Because the analysis and interpretation of drought frequency curves fol- 
lows that heretofore described for flood frequencies, two illustrations will 
suffice. Fig. 10 shows a group of drought frequency curves for the Milwaukee 


River at Milwaukee(13) with durations varying from one day to half a year. 


Examples of the interpretation of the curves can be seen from the two circled 
points. These indicate that the mean daily discharge will fall below 27 cfs at 
average intervals of two years, while the average weekly discharge will fall 
below 27 cfs only at average intervals of five years. Another method of rep- 
resenting this type of data is given in Fig. 11. Here minimum flows of the 
West Branch Mahoning River in Ohio(11) are plotted against duration. Read- 
ing the graph, it can be seen that the minimum mean monthly discharge of 
record equaled 0.05 cfs/mi2, but at the same time 75 percent of the years had 
a mean monthly discharge of at least 0.07 cfs/ mi2 and 50 percent had 0.095 


cfs/mi2, 


Frequency Analysis of All Flows 


Rather than developing frequency curves from only small segments of a 
streamflow record, analysis of the entire period of record has several engi- 
neering applications. Here, all of the mean daily flows constitute a popula- 
tion. Data are usually plotted on coordinates of discharge versus percent of 
time (of the period of record) during which any given flow rate was equaled 
or exceeded to form a flow-duration curve. 
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The flow-duration curve is familiar to many hydraulic engineers. Be- 
cause it encompasses the highest flow to the lowest, it is widely employed for 
planning and design purposes. Studies for surface water supplies, hydroelec- 
tric power developments, dilution and disposal of sewage and industrial 
wastes, construction diversions, and sediment transportation are based on 
flow-duration curves. Then too,the hydrologic characteristics of two drain- 
age basins can be compared by means of the flow-duration curves of each. In 
certain situations, it is even possible to treat the upper and lower portions of 
a flow-duration curve as flood and drought frequency curves. 

To prepare a flow-duration curve based on mean daily discharges, the 
number of discharges falling within predetermined intervals are counted. 
Discharges are divided by the mean discharge forming a dimensionless dis- 
charge ratio which facilitates comparison with other streams. The ratios 
are then plotted on a logarithmic ordinate and the cumulative time percent- 
ages on a probability abscissa. The flow-duration curve is sketched between 
adjacent plotted points. Log-probability paper enables a wide range of dis- 
charges to be plotted on one sheet and at the same time reduces the curvature 
of the flow-duration curve. 

Other than daily discharges can be used in developing a flow-duration 
curve; however, the shorter the period selected, the more definitive the 
curve becomes of the streamflow variations. Daily, monthly, and annual 
flow-duration curves are shown in Fig. 12 for an Ohio stream.(10) Notice the 
greater variability of flow displayed by the curves based on the shorter time 
periods. 

Curves based on mean daily discharges of four California streams are 
presented as Fig. 13. It can be seen that marked differences in flow distribu- 
tion exist among these streams, particularly regarding the highest and lowest 
flows. An investigation by Todd(27) of these and other California streams re- 
vealed that climatic and basin differences were strongly reflected in the 
shape of flow-duration curves. In mountain basins, where high flows result 
primarily from spring snowmelt runoff, the longer duration of these flows 
forms a relatively flat slope (Salmon River, Fig. 13). On the other hand ba- 
sins experiencing high flows from large short-duration rainfall show a rela- 
tively steep upper slope (Huasna River, Fig. 13). In a similar manner the 
shape of the lower portion of the curve is indicative of the ground-water stor- 
age capacity and aquifer permeability.(9) Curves for Arroyo Seco, Huasna 
River, and Thomes Creek in Fig. 13 show a pronounced steepness in the lower 
portions which is typical of basins containing only limited ground-water stor- 
age. In contrast, the flatter slope of Salmon River indicates a large, stable 
base flow. 

One limitation of flow-duration curves is that small changes in slope dis- 
guise important variations in discharge distribution. Because a flow-duration 
curve is the integral of a frequency distribution, slope variations of a flow- 
duration curve appear as variations in ordinates on a frequency distribution. 
Three such distributions are shown in Figs. 14, 15, and 16. Fig. 16 is based 
upon the same Salmon River data given in Fig. 13 and points up the greater 
sensitivity of the frequency distribution. 

Fig. 14 illustrates the common skewed distribution encountered in hydro- 
logic data. It should be noted that logarithms of discharge appear on the 
abscissa; if discharges were plotted on a linear scale, the skewness would be 
more strongly emphasized. The dashed curve is a theoretical Pearson Type 
III frequency curve(14) closely fitting the observed distribution. For 
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| comparison Fig. 15 shows the frequency distribution of a stream located not 
far from that of Fig. 14. The previous curve fitting is not applicable because 
of the different distribution of low flows. This curve does not approach zero 
frequency at the lower end; instead the low flows form a secondary peak lead- 
| ing to periods of no flow, which create a discontinuity on the logarithmic dis- 
| charge scale. The Salmon River distribution in Fig. 16 with its double maxi- 
mum displays a third type of frequency distribution. From a study of the 
seasonal flow distribution, it was found that the lower maximum was caused 
by the sustained base flow of uniform magnitude occurring during the fall and 
that the upper maximum was associated with runoff from winter rains sup- 
plemented by spring snowmelt runoff. 


CONCLUSION 


The purpose of this paper is to promote an understanding of frequency 
analysis of streamflow data. Emphasis was placed on concepts and advantages 
of such procedures; illustrations were inserted wherever possible to clarify 
points. References provide ready access to the omitted computational 
procedures. 

Whenever an engineering problem calls for information based on quanti- 
ties of streamflow, whether they be high, low, or all flows, frequency analysis 
can usually be helpful. It is hoped that a greater awareness of the analytic 
procedures, results, and limitations will continue the trend of statistical an- 
alysis becoming a more valuable tool for the civil engineer. 
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BUTTERFLY VALVE FLOW CHARACTERISTICS 
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J. C. Williams, Jr.,3 Junior Members, ASCE 
(Proc. Paper 1167) 


ABSTRACT 


Flow coefficients are given for free and submerged discharge as well as 
two cases with enclosed flow piping. The effects of type of installation, blade 
Shape and closure angle for control-type valves (using water) are demon- 
strated. Correlation of analytical relations with performance provides a basis 
for predicting incipient cavitation. 


INTRODUCTION 


Butterfly valves installed in pipes or penstocks are normally used for one 
of two basic services: either as closure or “on-off” valves, or for flow rate 
control or “governing.” Valuable information for fully-open service has been 
presented elsewhere.(1)4 Nearly all of the present discussion is limited to 
the flow of water through control-type valves for near-closed to near-open 
positions. 

The three major geometric variables affecting flow performance, in order 
of relative importance, are: (a) piping arrangement, (b) shape of blade, and 
(c) closing angle. For most applications flow characteristics are a function 
of geometry only. Available information on flow coefficients is limited and 
specialized. The first attempt at collecting and collating published data was 
made by Cohn,(2) but major variables were not delineated. 

This paper includes a demonstration of the effects of geometric variables 
and presents a basis for the analytical determination of flow coefficients. 
Also presented is a basis for the prediction of anticipatory minimum pres- 
sures related to cavitation inception. 

Note: Discussion open until July 1, 1957. Paper 1167 is part of the copyrighted Journal 
of the Hydraulics Division of the American Society of Civil Engineers, Vol. 83, No. 
HY 1, February, 1957. 

1. Associate Prof. of Civ. Eng., Lehigh Univ., Bethlehem, Pa. 

2. Asst. Prof. of Civ. Eng., Univ. of Washington, Seattle, Wash. 

3. Senior Project Engr., CDC Control Services, Inc., Hatboro, Pa. 

1, Numbers in parentheses refer to references listed at the end of the paper. 
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Definition of Terms 


Installations with equal internal diameters of valve body and upstream 
pipe have been investigated, with all variations in piping arrangement made 
beyond the valve. In figure 1, Station @) is sufficiently upstream from the 
valve as to insure unaffected, parallel streamlines (uniform flow). Condi- 
tions at Station are variable. Station is generally within the valve 
body and represents the cross-section of the jets where uniform flow takes 
place. For the sake of correlation, the flow formula adopted by Cohn is used 
throughout: 


Q=C,D* Vg vaH (1) 


where Q is the discharge in cu. ft./sec., 
Cg _ isa coefficient of flow, 
D is the internal diameter of the valve body in feet, 
g is 32.2 ft./sec./sec., 
AH is the net difference in total head between Station @) and G) ; 
in feet of the fluid flowing; i.e., the valve loss. 


In order that the reader can immediately appreciate the effect of piping 
arrangement on the flow coefficient, the following flow coefficient curves are 
presented in Figure 2 for a specific valve. Station @) for each of the curves 
respectively is: 


a) Downstream from a 4''x 6''diffusor, 9''long; 

b) in a pipe of same diameter as Station and the valve, 

c) in the atmosphere beyond live jet (water discharging into air); and 
d) in still water beyond jet (water discharging into water). 


Referring again to Figure 1, “left” applies to the downstream half of the 
blade and “right” to the half projecting upstream of the hub. For each value 
of the angle of inclination of the blade, &, (0° at wide open position), there is 
a corresponding specific value of x and y, left and right. The % and y dimen- 
sions may be affected by viscosity and might therefore vary with both ao and 
a Reynolds number. 

This paper is restricted to a consideration of incompressible turbulent 
flow. 


Basic Analysis - Free Discharge 


Flow between Stations @) and 2) is accompanied by an increase in velocity 
with distance, a highly efficient process. For all installations reported here, 
the head or energy loss between these two stations due to the flow accelera- 
tion is therefore neglected as being relatively minor. 

In true free discharge the pressure at Station @) is atmospheric. Since 
uniform flow occurs at this station (so long as Zy, does not project beyond the 
valve body), the pressure distribution in the jets at this station is hydrostatic 
and at equal elevation the velocities of the two jets are equal. 

If the areas of the two jets at Station (2) are known, it is then possible to 
estimate the jet velocity head (in this case, the head or energy lost from the 
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System) for a given blade setting, without recourse to an actual flow 
measurement. 

The jet area at Station Q) is bounded by a circle and an ellipse. The jet 
area on one Side is therefore equal to mDy/4. Since the jet area is a function 
of the variable y, measurements of the jet thickness in a two-dimensional 
flow model, representing a narrow section of the center part of the blade, 
should provide a basis to understand and predict prototype performance. 


ee 


Two-Dimensional Flow Distribution Study and Prototype Verification - Free 
Discharge - 6-Inch CDC Valve 


A two-dimensional plastic flow section (Figure 3) was constructed (cross- 
section 8"'x 1") and measurements of yy, and yg were made for the four blade 
shapes shown in Figure 4. 

The length from the centerline of the blade to the end of the valve body is 
5D/8 for both the 4-inch and 6-inch CDC prototype valves (for consistency 
this same length was used for all four blade shapes in the two-dimensional 
free discharge tests). 

Since the pressure at both left and right jets at Station @) is zero, both jet 
velocities are equal. It follows then that 


Var= 1. / Ye (2-2) 
and Voee e/a (2-b) 
then a/, 2 Ya, (2-c) 


where qj, + q@R = QT, total flow in cu. ft./sec./unit width. For three- 
dimensional flow, the area of one jet is mDy/4, rather than y as above, but the 
discharge ratio is the same - 


nay 
Q, /Q;= Y.+Yp 


where + = Q7, total flow in cu. ft./sec. 

Pays eee the first to make an issue of the fact that the quan- 
tity of flow leaving each side of the blade is not equal. 

To substantiate Gaden’s contention, free discharge tests were run on the 
two-dimensional apparatus and measurements of y and q were made at vari- 
ous blade angles, using a valve body length equivalent to the prototype valve, 
the 6-inch CDC valve. Measurements of flow distribution were made also on 
the prototype valve. The results appear in Table I, and were of sufficient 
precision and promise to encourage further study. By 

In the prototype tests, a limiting minimum Reynolds number (NR = V1 in ; 
where Vj is the average approach pipe velocity, D is the u stream pipe an 
valve diameter and p is the fluid kinematic viscosity in ft.“/sec.) was foun 
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FIGURE 3 


TWO DIMENSIONAL JET APPARATUS 
SHOWN SET UP FOR THE \" CDC VALVE FOLLOWED BY A DIFFUSC 


(6" Free Discharge Piece and Extension for 
"Continuous Pipe" are shown at bottom of photo). 
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for a constant flow. coefficient, CQ, as depicted in Figure 5. The prototype 
flow ratios of Table I represent values of Np in excess of these minimums. 
Referring again to Figure 5, the ratio Q1,/QR ceased to be constant once the 
NR fell below the minimum for a constant CQ. Figure 5 should be a guide for 
free-discharge conditions in general, but the effect of a different blade shape 
and/or valve body length would be unknown factors with another valve. 


Free Discharge and Submerged Flow Characteristics 


For free discharge, then, there is a definite minimum Reynolds number, 
Np, for a constant flow coefficient. This is due in part to the fact that y, and 
particularly 4, (Fig. 1), may be influenced by boundary-layer effects at 
lower Reynolds numbers. 4;, might thus be contained within, project beyond 
or coincide with the termination of the valve body at different NR values. 
Also, the jet may be affected by distortion of the spray in the lateral expan- 
sion beyond the uniform flow at Station . Therefore, if the ratio of jet 
thickness is upset, both Cg and Qy,/QrT will be affected. 

As stated previously, the head loss from Station to @) due to accelera- 
tion should be insignificant inasmuch as an increase in velocity obtains. 
With h being the piezometric head at a given station (pressure head plus ele- 
vation), the equation for free discharge into the atmosphere with hg = 0 be- 
comes: 


V Vie 
agth ag 4H 


2 
(3) 


With submerged flow (flow of water from valve discharging directly into a 
large body of water), the equation is assumed to be: 


2 


a5 h-h,=aH 4 


where hg is the depth of still water above the valve centerline. Free dis- 
charge tests were run on the 6''CDC valve, and free discharge and submerget 
tests were run on the 4" CDC valve. The 4" valve test results, given in 

Table II and plotted in Fig. 2, show that there is a distinct difference between 
the free discharge and submerged flow values. This difference had been in- 
ferred in studying available data, listed in Table III and plotted in Fig. 6. 

The Schmidt and Davis data (air flowing into air) are given by Cohn, in he: 
Fig. 7, p. 883, but details are lacking. The Genissiat data(4) are for a valve 
body which converges to 84% of the upstream pipe diameter, between the cen- 
terline and the end of the valve body, and are included for comparison only. 
The Dow data are given by Gaden(5) and apparently are for the 78-inch valve 
of Fig. 8, p. 667 by Dow.(6) (The original Genissiat data was given by 
Laurent).(7) 

Note that the data for the 4-inch valve tested under submerged discharge 
conditions corresponds quite closely with the data of Schmidt and Davis, ex- 
cept in the near-closed positions. This discrepancy can be explained by the 
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TABLE I 


6-INCH CDC VALVE DATA - FREE DISCHARGE 


____. _Two-Dim, App. Measurements ____ Prototype Valve Data 
x Sa ae eh "Pipe. : (meas) cnet 
pe oe 4 Tyee ly 2 NR ara (Fig. 5) 

oO ~—-.498 500 5.6 x 10? 3 4.8 x 105 

5° 50k, 250, = 6.1 - 4.6 
BOp ews 560 2569-546 - 4.3 
15% 568 -561 = 55.0 - 4.2 
20> S70 2568 349 - 348 
25° 555 2562 (3.6 - 3.6 
30° SS 255303 0.560 303 
35° - - - 0.550 3.0 
4o° = 6556 055 | Os 0.510 257 
ys? ho 5uB 046 0.535 2.3 
50° 4537 534 08 0.532 1.9 
55° 513 517 049 0.528 vee 
60° —_.520 50h 067 0.525 1.0 
65° = 53 52l (Oe 0.523 - 


“at sufficiently high Np to obtain constant Q-ratio; required minimum Np clos 
that per constant Cg of Fig. 5. These are averages of flow measurements on e. 
side of the blade, i.e., Cy vs Qp and Qr vs Qpe 
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TABLE Iil 


FLOW COEFFICIENT Cg 


Free Discharge Submer zed Free Discharge 
C= Ces: Schmidt Davis Génissia 
a Lng 6" g (Converging) Dow 
200 mm 
0° - - 1.30 - 0.7h 0.67 
oe - 0.63 = - - = 
10° 0.71 0.620 0.9 0.90 0.69 0.66 
oh 0.635 0657 - - 0.61 0.61 
20° 0.565 0.528 0.65 0.64 0.53 0.55 
AG? 0.190 0.476 - - - 0.47 
30° 0.420 0.405 0.45 Ohh 0.38 0.40 
35° 0.355 0.35 - - - 0.33 
hor 0.285 0.285 0.31 0.29 0.26 0.27 
ie. 0.223 0.232 - é = = 
50° 0.166 0.183 0.21 0.19 0.15 - 
55° 0.120 0.13h * 2 2 4 
60° 0.083 0.088 0.11 0.12 0.09 - 
65° 0.053, + = “ : s 
70° - - 0.06 0.07 = a 
75 x - = - - = 
80° x x 0.0 0.02 e = 
85° x x x x x Z 
90 x % x x x = 


Closure 75° Ts 90°? 90°? 90° U 
Upstream ¢ 6 he? unknowm unknown 0283? n 
tmax/D 0.269 0.225 unknom  unlnomm 0.225 k 
Fluid water water Air Air water n@ 
0 
w 
n 
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large difference in closure angle. Good correlation was also obtained be- 
tween data from the 6" and 4" free discharge tests. These data, fortified by 
the other free discharge data, indicate a separate trend of their own. 

The free-discharge and submerged characteristics appear to be related, 
but are distinctly different at the near-open positions. In submerged flow, 
the fluid fills or nearly fills the void between the two jets up to the blade, re 
sulting in an expansion of the given fluid within the same fluid. In free dis- 
charge of a liquid, the void is normally filled with air. For this reason, it 
would be logical to assume that the flow characteristics might not be identi- 
cal even though the installations are quite similar. This assumption is born 
out to some extent by the fact that the free discharge analysis was substanti- 
ated more closely by prototype data than was the submerged case. 


Continuous Pipe Flow Characteristics 


By a continuous pipe is meant a constant diameter D between Stations @) 
and . With the 6"'CDC valve it was found that for this type of installation 
hydraulic grade (or piezometric head) lines unaffected by the valve could be 
obtained as close as 2D upstream and 5D downstream from the valve body. 
In the near-closed positions the distance requirements are less. In general. 
a distance of 2-1/2 D upstream can be recommended as safe. Lack of suf- 
ficient data prohibits any recommendation concerning the safe downstream 
distance. 

The continuous pipe tests were conducted at a wide range of Np, in many 
instances well below the minimums set forth in Fig. 5 for free discharge. 
The lowest pipe NR reached was 8 x 104 (ato 4 = 45° with the 4" valve, and < 
Oo = 65° with the 6" valve). The flow coefficient for o 9 > 30° was relatively 
constant and no obvious viscous effects were noted. This is entirely logical 
since the head loss in the expansion beyond the valve is practically independ 
ent of boundary friction and due primarily to the flow expansion from Statior 

to (3). 

For this piping arrangement, since Vj = V3, 


Rehyee cl (5) 


The data of Table IV, plotted in Fig. 7 are in close agreement for a9 > 2 
even though the blade shapes and closure angles differ. Flow coefficients fc 
the near-open positions are influenced by the relative blade thickness and th 
blade shape. According to Netsch and Schulz,(8) viscous effects are noted 
with Oo = 0° to 20°. Attainment of a high degree of test accuracy within thi: 
range is difficult because of the relatively large discharge rates and small 
head losses involved. 

In all data presented in Table IV a deduction from the overall loss has 
been made for equivalent pipe friction to and from the valve body for purpos 
of correlation. 

The Netsch & Schulz data(9) are all for an Np of 3.5 x 109. 

The Dickey and Coplen data(10) for Disk “I” and Disk “II,” with D about 
13.5 inches, are for air velocities from 500 to 5000 ft./min. A great variet 
of downstream pressures were used. They state that Co is unaffected by si 
or velocity. Disk “I” was a thin, flat plate with short hubs only at the outer 
edges. Disk “II” was of uniform thickness with edges rounded to 6% of D. 
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| TABLE IV 


FLOW COEFFICIENT, Cg 
CONTINUOUS PIPE DATA 


Netsch & Dickey- Dickey- Genissiat Square 
cDne CDC Schulz Coplen Coplen (Converging) Damper 
3 yu" pg 6" g 10 cm DiskI Disk II Dickey- 
Xe Blade "Bp" (Thin) Coplen 
0° 1.65 - 5.97 = = = = 
> = oe - - - ~ 3.54 
10° 1.38 1.23-1.35 1-78 yey 1.56 2eh7 2.35 
15° ~ - - - - 1.59 S 
20° 0.92 0.89-0.99 0.01 1600.5 1503 1.10 1.29 
250 - = = = 2 5 s 
30° 00592 0.56) 0.582 0663 0.625 0.62 0.767 
35° - - 0.440 s = = x 
40° 0.354 0.330 0.351 0239 0638 0.38 OLS 
45° 0.272 0.260 0.269 - ~ ~ ~ 
50° 0.208 0.19) 0.202 0.22 0.22 0.21 0.272 
55° ad a 0.11 ™ or = a 
-60° 0.106 0.092 0.100 0.120 0.125 OL 0.147 
65° 0.064 0.059 0.060 - - - - 
70° - - 0.020 0.062 0.063 - 0.0685 
5° x = = im = * 
80° x = x 0.028 0.037 S 02027 
g5° ibs et 5 © 2 = : 
90° x x x x x x x 
a en, ol RO 
Closure 75° -77.5° 722 90° 90° 90° 90 
Upstream 9 2 2 2.5 unknown unknom ~~1l unknown 
D'nstream B 5 5 3.5 unknom unmom ~7 poe ak 
+ 00269 04225 0.12 0.01 0.17 06225 ° 


Fluid water water water Air Air water Air 
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Their tests were run at 80°F. (The square damper had the same profile as 
Disk “I”, and is included for comparison only). 

The Genissiat data is for a D of 155 mm, converging to 84% of D from the 
centerline of the valve to the end of the valve body. These data are included 
for purposes of comparison. 

Lack of consistency and/or necessary details precluded use of data for 
Gaden’s Disk “A” or for data presented by Keller and Salzmann.(11) 


The 4" CDC Valve Followed by a Specific Diffusor 


In this series of tests the 4''CDC valve was preceded by 4''diameter pipe 
and followed immediately by a standard 9"' long diffusor and a 6"' diameter 
pipe. The basic data are listed in Table V and plotted in Fig. 2. The flow 
coefficients for this arrangement are nearly equal to those for a continuous 
pipe. It should be stressed, however, that this comparison should not be ex- 
pected for any other combination of valve and diffusor. 

As with the continuous pipe tests, no viscous effects were noted for Np as 
low as 8 x 104 (at 45° and 65°). A deduction from the overall loss has been 
made for pipe friction, but no deductions were made for the wall friction 
losses through the valve body and the diffusor. 

For this arrangement: 


(6) 


Sufficient piezometer stations were established to insure adequate deter- 
mination of gradients: 


Number of Number of 
Valve Installation Upstream Stations Downstream Stations 
ORECDC Continuous Pipe 5 8 
4"'" CDC Continuous Pipe 4 4 
4"" CDC Diffusor 4 10 


General Analysis for Enclosed Flow 


As argued previously, the loss of head or energy between Stations (1) and 

in addition to pipe friction should be relatively minor. In contrast, there 
is an abrupt decrease in velocity with distance between Stations @) and (3). 
Separation occurs within the fluid for enclosed flow conditions (such as with 2 
continuous pipe or a diffusor, or submerged). 

The high intensity of shear in the vicinity of the separation surface 
and the eddy rollers formed between the live jets are manifested in a 
non-recoverable head or energy loss. The initial point of expansion is at 
Station (2), for both jets, and is not far from the initial point of separa- 
tion at the edge of the blade tip. In general, it has been observed that - 
the piezometric head within the zone of separation (here, between the 
live jets) is constant, and equal to the piezometric head at the point of 
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TABLE V 


4=INCH CDC VALVE FOLLOWED BY A DIFFUSOR 


Notes: 

lk" $ CDC valve followed immediately 
by a 9", 4" x 6" diffusor (internal angle, 
12.5°). Uniform hydraulic gradient for data 
obtained extends to about 8 valve diameters 
downstream as opposed to only about 5 valve 
diameters for continuous piping; two valve 


diameters were used for the upstream station. 
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separation.(12) Within the live jets at Station (2) the pressure distribution 
should be hydrostatic (constant piezometric head), since uniform flow occur 
at these sections. It is assumed that the piezometric head is the same in th 
live jets at the separation surface at Station (2) and within the adjacent zone 
of separation. With this assumption, it would be within the same contracting 
surface of separation. With this assumption, it would be reasonable to con- 
clude that there is some path between the left jet and the right jet at Station 

over which the pressure distribution is hydrostatic. To illustrate, if the 
longitudinal section of Fig. 1 is placed in a horizontal plane (for liquid flow), 
the pressure within the left jet at Station @) is thus equal to the pressure 
along some connecting path, through the zone of separation to and within the 
right jet at Station (2), in this instance the pressure within the two jets woule 
be constant and equal. 

As was the situationfor free discharge, it then follows that Var, equals 


V2 Determination of the velocity distribution at the point of separation 


Rr 
would not be practicable. The relative magnitude of the velocity at Station ( 
however, can be expressed in terms of the uniform flow width, y. Inasmuch 
as the surface of separation originating from a fixed point (here the blade ti 
is usually fixed in space for all but low values of Np, y should be constant fc 
a given blade setting. 

The continuous pipe and special diffusor installations were duplicated wit 
the two-dimensional test device and values of y were measured. In these 
tests, the magnitude of y was not noticeably changed by varying the discharg 
Differentiation of the live jets was facilitated by injecting dye or air bubbles 
downstream from the blade. 

The flow expansion between Stations @) and GB) is extremely abrupt and t 
effect of boundary-layer development, or the boundary shear stress involvec 
is comparatively negligible. Once the jet thickness is known, it is possible | 
appraise the energy loss due to the expansion alone, by means of the princip 
of conservation of momentum. It is assumed throughout that the valve loss" 
very nearly equal to the loss due to flow expansion, with the loss along the 
boundary between Stations (1) and (3) and within the fluid between Stations (1 
and Q) being relatively minor. 


Comparison of Actual and Analytical Flow Coefficients 


The jet area has been shown to be a function of y (see Fig. 1). Also, 
Vi= Vr for all piping arrangements, so long as any boundary alignment 
change occurs downstream from the distance 2. 


Then QL/QR = YL/YR; 
or Qf = V2 x (the area of the two jets). 
Or, QT = (V2) (tDya/2) 
and Ay/Aj = D/2ya, 

where yg is the average of yz, and yp, 

Aj is equal to Q/V2, 
and Ay is the area of the valve body. 
Note that Ay is identical with the area at Station (4) for all cases studied. 
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It is now possible to write an equation for the conservation of momentum 


tween Stations (2) and (3), with the basic relationship established for the 
low areas. 


For free discharge, since AH = V22/2¢ = V12/2¢ + hy, it may be shown 
hat 


rh2V2) (7-2) 
AV/A, 


For a continuous pipe, theoretically AH = (V2 - V3)2/2¢ (the Borda-Carnot 
loss), from which 


Cg= 


+/(2/2) as) 


For the particular 2 = 3 ratio diffusor tested, assuming that 
AH = (V2 = V3)2/2g, 


(7-c) 


Both the free discharge and the continuous pipe data (particularly the lat- 
ter) provided good agreement with these equations for Cq using values of Ya 
obtained from the two-dimensional apparatus. The diffusor data was consist- 
ent, and when arbitrarily corrected as follows, provided equally good 
agreement: 


_whev2) (1-ce) 


oA 4 
b2ha. 9 


By coincidence, the area ratio Ay/Aj is equal to D/2yg for both the two- 
jimensional and three-dimensional cases. Although the two-dimensional 
yalues for Ay/A; using D/2ya were fairly close to those computed from mea- 
sured values of Cq and Eq. (7-a), (7-b), and (7-cc), quite close correlation 
‘or all installations was obtained by adjusting the two-dimensional values as 
ollows: 


5 (8) 


A\/Ay= La(D/2y 


| 
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Eq. (8) applied equally well to both the 4-in. and 6-in. CDC valves. Inu 
ing D/2y, as the area ratio, it had been assumed that the flow was strictly 
rectilinear at Station Warpage of flow filaments was evident in the 4-i 
CDC valve free discharge tests, particularly on the “left” side. The result 
this instance was a much wider spray pattern in the vicinity of the “left” ar 
“right” side than in the larger periphery between. A picture illustrating b 
the warpage and the uneven spray is presented in Figure 8. In this arrangt 
ment, a special test not otherwise discussed in this paper, the “left” side c 
the blade is below the horizontal shaft. 

The 4-inch CDC valve area ratios are given in Table VI and those for t 
6-inch CDC valve in Table VII There is good agreement between all 
columns, showing the validity of the analytical equations. Note that the are 
ratios are little influenced by downstream conditions. There is not much ¢ 
ference between the 4-inch and 6-inch CDC valves, since the effects of blat 
shape and closure angle on Cq are obvious only at the near-closed or near 
open positions. 

To illustrate the possible use of the two-dimensional apparatus in predi 
ing performance, some values of Cg for continuous piping and free dischar 
for the 4-inch CDC valve were computed from apparatus data using Eq. (7- 
(7-b), and (8), and are plotted in Figure 2. 

Reference data is given for continuous pipe arrangements in Table VIII. 
The two-dimensional “Synthetic Blade” was very similar in shape to the 
blade used by Netsch and Schulz. A special flat blade was fabricated for tt 
two-dimensional studies comparable to the Dickey-Coplen Disk “I.” (See 
Fig. 4). 

In passing, it may be noted in Fig. 1 that two separation surfaces are 
shown for the “right” jet. For a blade with a relatively wide tip, the point 
separation will be at one or the other edge, depending upon the relative bla 
angle. The sudden, slight shift of CQ trends for some valves is due to this 
transfer of the separation point. 


Prediction of Incipient Cavitation 


With Vj, hy, hg, and oo given for an installation, values of V9 and henc: 
hg can be calculated once the magnitude of Ay/Aj is known. By using C 
values from flow tests in equations such as (7-a) and (7-b), it is believed t 
a reasonable approximation of Ay/Aj should be obtained from which hg car 
computed. An anticipated hg for specific conditions is outlined in the sam) 
calculations of Table Ix. 


Fully Open Position 


In the fully open position (a9 = 0°) the point of separation has been foun 
by Mueller(13) to depend upon the blade shape and the Reynolds number in. 
volved. The variation in Cg for the 6" CDC valve, Table IV for 10° and 20 
might have been the result of a shifting point of separation. As stated pre 


viously, Netsch and Schulz (8) noted viscous effects for continuous piping w 
Qo = 0° to 20°, 
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TABLE VI 


COMPARTSON OF AREA RATIOS FOR )" CDC VALVE 


as computed by using analytical equations for C,, and as 


obtained from two-dimensional flow test measurements. 


—_ From Flow Tests Two-Dim. Measurements 
ae ar i fa. (8) 
a 
a Eq.(7-a)* Eq.(7-a). Eq.(7b). Eq.(7-cc) : 
Submerged Free Cont. Diffusor Free Cont. Diffusor 
Xx. Discharge Pipe Discharge Pipe 
0° = = = = = - abaeen 
= = 2 x = 2 % x 
10° 1.31 1.56 1.80 ~ ~ ~ 1.84 
5° 1.50 1.75 - - 1.96 1.91 2.02 
20° ee 1.97 raya lan arse le 2,18 2.18 2.2 
25° 2.14 2.26 2.50 - 2.51 2.50 2.58 
30° 2.62 2.64 2.89 2.94 2.944 2.92 3.02 
35° Beet 3.12 3042 - 3.53 3.50 3.65 
40° 4.19 3.89 4.18 = 39 4.30 he2h held 
45° 5.70 4.98 5615 508 5.37 5.28 5.58 
50° 7.9 6.68 6.50 6.90 6.72 6.55 6.50 
55° 211 9.25 8.41 - 8.72 8.40 8.60 
60° 15.9 13.4 11.58 12.65 ie; 9 dl, eee lee! 
65° etee 20.9 17.35 20.25 19.4 18 .2 Zine 
70° 31.8 - - - - - ~ 


ee ea ees ea ie ETE EL 
3 Assuming Eq. (7-3) for Free Discharge can be applied to Submerged Flow. 
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TABLE VII 


COMPARISON OF AREA RATIOS FOR 6" CDC VALVE 
as obtained by using analytical equations for Co» and 
as obtained from 


two-dimensional flow test measurements 


From Flow Tests Two-Din. Measurements 
3 0.85. 
A/a; A/h; = (Be) 3 Ea. (8 
rent a Free Cont. 
Ae Discharge Pipe Discharge Pipe 
10° 1.80 - 1.65 a 
157 1.95 = 1.90 1.97 
20° 2.11 - 2.15 2.22 
25° 2.34 es 2.51 2.58 
30° 2.7h 2.97 2.90 2.98 
35° 3.23 3.53 3.40 3.57 
40° 4.79 5.28 5.05 5230 
50° 6.07 6.73 6.35 6.60 
55° 8.30 849k 8.10 Beh5 
60° 12.55 3.1 12.0 11.60 


65° - 19.9 2125 18.2 


SCE 
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TABLE VIII 
COMPARISON OF AREA RATIOS FOR VARIOUS VALVES 
as obtained by using analytical equations for Cg, and as 


obtained from two-dimensional flow test measurements 


All for Continuous Pipe 


From Flow Data Two-Dim. Measurements 

A,/A;, from Eq. (7-b) Ay/A, = 162 ae Bq. (8) 
Netsch Dickey=- Square Synthetic Flat 

& Coplen Damper Blade Blade 
Schulz DiskI_ De ~Ce (le. & Se) (Disk 1) 

- - 1.4 - - 

- 1.60 1.61 - 1.58 

- - - 1290 1.84 
2.10 2.11 2.10 2.16 2.14 

= = - 247 2047 
2.91 2.7% 2.85 2.87 2.90 
3.53 = - 3.35 3-7 
4.17 3.85 4.12 4.03 e10 
Seth - ~ he97 5.00 
6.50 6.05 6.20 6.20 6.20 
8.90 - - 8.20 8.00 
12.10 10.3 10.65 11.8 10.5 
19.55 - - 21.5 1.9 

- 19.0 21.65 - . 
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TABLE IX 


SAMPLE CALCUIATIONS — ANTICIPATED hy 


The h-inch CDC valve installed in horizontal “continuous pipe" with 
hz = 10 fte, Q = 2.00 cfs and o, = 35°, compute hy. 


From Fige 2 Cg = Ooh6 

From Eqe (1) 0H = O°/(Cg?D4g) 
AH = (2.00)?/(0446)7(06333)4(32.2) 
AH = 7.7 ft. 

Therefore hy = AH?r h3 = 57e7 fte 


Using Eqe (7-b), Ay/Ay = 7 +1 = 3.2 
a ama ot 


With By # Ay = SS 
Then 4, = 0 lad = 0.0255 sq. fte 


V; = 2.00/0.0255 = 78. ft./sece 
v,*/2e = 95. ft. 


V, = _2.00_ = 23.0 ft/sec. 


TW 
v,°/28 = 8.2 ft. 
Since r°/ 2g th, = ei 2g + hoy 
ho = 862 t 5767 — Mel) = — 29.5 fte 
For 60°F water service, vaporization is probable and cavitation effe 


Note that hy and hy are related directly to h3 (although hy is limit 
vaporization of the liquid) as shown in the tabulation below: 


h3 ft. hy ft. hy ft. 
10 57e7 -29.5 
20 6707 19.5 
30 T71e7 - 965 


ko 8767 +0.5 
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SUMMARY 


The effect of piping arrangement on the flow characteristics of a specific 
lve has been demonstrated. Evidence has been submitted to prove that 

se, submerged and continuous piping discharge arrangements have separate 
d distinct flow characteristics. For the valves discussed in this paper, the 
ww characteristics for a given type of installation were not appreciably in- 
‘enced by either the shape of the blade (except at the near-open positions) 
by the closing angle (except at the near-closed positions). 

Analytical relationships have been presented. Measurements from a two- 
mensional model have been combined with those for prototype valves to 
monstrate the basic validity of these relationships. As a by-product of this 
idy, it may be noted that relatively good anticipatory values of flow coeffi- 
ents for other types of installations could be obtained by means of inexpen- 
ve two-dimensional tests using water, air, smoke, a polarized or diffracted 
lid, etc. 

Probably the most important use of the relationships offered will be in 
lculating anticipated conditions for incipient cavitation. 

It is felt that the accuracy attained in the studies reported is adequate for 
rmal application of the results. 
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Discussion of 
“TIDAL CURRENTS AT INLETS IN THE UNITED STATES” 


by Joseph M. Caldwell 
(Proc. Paper 176) 


JOSEPH M. CALDWELL, ! M. ASCE.—Mr. Baines has presented a rather 
nteresting analysis of one facet of tidal current action at inlets. It was, of 
Ourse, the writer’s intention in the paper to present a very simple hydraulic 
lassification of tidal inlet based on a minimum of data. In this sense, Mr. 
3aines’ analysis goes beyond the simplified concepts presented by the writer. 

Mr. Baines prefers to treat the writer’s class 2 and class 3 inlets as 
xamples of a general type of inlet. The class 3 inlet (the inlet with an inade- 
uate entrance) is not influenced by the length of the contiguous bay or estu- 
ry. On the other hand, the class 2 inlet (short inlet with adequate entrance) 
Ss so definitely related to the length of the contiguous bay or estuary that, in 
he writer’s thought, the distinction between the class 2 and class 3 inlets 
hould be preserved. Methods of treating two or all of the classes 1, 2, and 

of inlets by a single mathematical concept (Einstein,(1) Dronkers)(2) in- 
ariably leads to complications; the writer has attempted to assist the user 
ho is not in a position to utilize the more complicated methods in approach- 
ng a problem dealing with tidal inlets. 


REFERENCES 
. Einstein, H. A., “Computations of Tides and Tidal Currents— United States 
Practice,” Proc. ASCE, Vol. 81, Separate No. 715, June 1955. 


. Dronkers, J. J. and J. C. Schonfeld, “Tidal Computations in Shallow 
Water,” Proc. ASCE, Vol. 81, Separate No. 714, June 1955. 


7 Chief, Research Div., Beach Erosion Board, Corps of Engrs., U. S. Dept. 
of the Army, Washington, D. C. 
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Discussion of 
“EXTENDING STREAMFLOW DATA” 


by W. B Langbein and C. H. Hardison 
(Proc. Paper 826) 


W. B LANGBEIN** and C, H. HARDISON,** A.M., ASCE.—The authors 
eply appreciate the very profitable and constructive discussion and criti- 
sm of their paper. We are particularly indebted, not only to Mr. R. D. 
20drich’s discussion, but also to many of his earlier works for numerous 
iggestions and findings of value. 

Messrs. Merriam and Schuleen introduce a number of items that unfor- 
nately could not be presented with the ample consideration they properly 
serve within limits of a discussion. They tempt us with all-too-brief pre- 
tations of new methodologies and ideas they have developed. We were 
ticularly interested in their use of punched cards to develop interstation 
lations. 

The authors are in substantial agreement with Mr. W. M. Snyder and 
knowledge his efforts to clear up several important details. 

Mr. Snyder mentions that “the use of the geometric mean to obtain a unique 
1e of relation between two variables has been mentioned by other authors.” 
e credited this device to Goodrich whose description of it was apparently 
out contemporaneous with that by Gumbel. 

The equations and discussion by Snyder on the variance of a forecast are 
eated in an earlier paper! where it was shown how the correlations between 
¢ flow of adjoining streams can be applied to the design of networks of base 
d secondary stream-gaging stations. 

Snyder recommends that we use the term “standard deviation of the re- 
duals” in lieu of the more usual “standard error of estimate.” Our usage 
confirmed by several texts, but to be strictly correct for small samples 
2 standard deviation of the residuals should be adjusted for the degrees of 
2edom lost in drawing the curve of relation. Thus on figure 1, assuming 6 
grees of freedom consumed by the nature of the graph shown, the standard 


ror of estimate should be 0.086 \ as or 0.0906 logyo units. The standard 
ror of a graph, Sg, is usually computed by the following equation: 


a ef ees (1 


ere s is the standard error of estimate computed as above, x-x is the 
viation from the mean of the x values, o is the standard deviation of the x 
lues, and n is the number of data used to define the relation. Eq. 1 pre- 
mes random data, whereas the data we have to deal with are time series. 


‘Hydr. Engrs., U. S. Geological Survey, Washington, D. C. 
, Langbein, W. B, Stream-gaging networks, Publication No. 38 de la 
Association Internationale d’Hydrologie, 1954. 
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This means our data have some properties that differ from the kinds of date 
usually treated in statistical literature. One of the more important of these 
differences is that our data are related to each other in sequence; that is, 0! 
day’s flow is somewhat related to that of the preceding and following day. 
The same is true of months, seasons, and years. For example, as observec 
by Johnson, deviations from a correlation tend to persist for periods of 
several months. Our record, therefore, consists of data more or less inter 
dependent; or, in different terms, our record represents fewer independent 
events than would be indicated by the number of items in the list. The actuz 
number of independent events, using monthly data, depends on the tendency 
toward serial correlation among the monthly data. Further research is nee 
ed on this point, but at this time we would judge that an estimate of the prop 
value of n in Eq. 1 could be taken as twice the number of changes in sign of 
the net deviations arranged chronologically (assuming plus and minus signs 
to be equal in number), but in no case to exceed the number of months. Thi 
rule applied to the Murder Creek example would yield 40 independent data i 
stead of 60. In the case of Rock Creek in Utah, the indicated number of inde 
pendent data is 36 out of a 72-month list. 

In usual practice, the index of correlation is computed from the ratio of 
the standard error of estimate (adjusted for degrees of freedom lost in defi 
ing the relation) to the standard deviation of the monthly discharges. How- 
ever, because of the tendency for similar months to cluster, it is believed 
that more realistic values of the correlation index can be computed by takin 
the ratio of the adjusted standard error of estimate to the standard deviatio 
of the monthly discharges about their respective monthly means. To be cor 
sistent, this standard deviation, obtained by pooling variances, should be ad 
justed for the 12 degrees of freedom lost in defining the monthly means. 

In the Murder Creek example, the coefficient of correlation of 0.93 show 
in figure 1 is computed from a standard deviation of the residuals (error of 
estimate) of 0.086 and a standard deviation of 0.231. The standard error ad 
justed for 6 lost degrees of freedom is 0.0906 as computed previously. The 
standard deviation of the monthly discharges, about the means of the respec 
tive calendar months has been computed to be 0.18. The standard deviation 
adjusted for the 12 degrees of freedom lost in using the means of the calen- 
dar months becomes 0.18\/ soar or 0.201. 

The index of correlation computed from these adjusted values is: 


2 
ie es nae 


0.201 


In this example the change in the index of correlation is not great, but in 
streams where the seasonal variations are more marked, as in the case of 
western snow-fed streams, it is of greater significance. In one case, for 
example, Rock Creek vs. South Fork Ogden River in Utah (70 miles distant 
the following were observed: 


Standard error of estimate = 0.12 logyg units 
Standard deviation of monthly discharges about the general mean = .41 


Pet oy (2) = 0.95 


41 
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However, it was also observed that standard deviation of the several 
months about the means of the respective calendar months is not much larger 
than the standard error of estimate. Computing the index of correlation rela- 
tive to the standard deviations (adjusted for 12 lost degrees of freedom) 
about the mean of the respective calendar months gives a non-significant 
value and a curve of relationship apparently offers no special advantage in 
this case. 

For a correlation between monthly flow of Rock Creek and that of Yellow- 
stone Creek, 18 miles apart, the respective factors are: 


Standard error of estimate (adjusted) = 0.050 logy9 units 
Standard deviation of monthly data about general 
mean = .41 
Standard deviation of monthly data about means of 
respective months (adjusted) = tally 
Index of correlation relative to general mean —.09 
Index of correlation relative to means of calendar 
months = OL 


Consequently, the flow of Rock Creek correlates satisfactorily with that of 
Yellowstone Creek but not with the flow of South Fork of Ogden River. 

Use of rainfall data: It is the practice to include rainfall data in extending 
stream-flow data wherever such use is beneficial. The authors suggested, as 
Mr. Johnson was kind enough to demonstrate, that adjustment of the correla- 
tion between the flows of two streams on the basis of the relative rainfall on 
the two drainage basins, might afford a superior arrangement. 

In extending a record in time there is always the question whether the re- 
lationship may shift due to such factors as a shift in rainfall distribution. In- 
cluding rainfall in the correlation or correlating with streams on either side 
normal to the prevailing storm tracks as suggested by Johnson can minimize 
the effects of shifts from this cause. The greater number of independent 
jata that are used, the more stable the correlation. 
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Discussion of 
“INTEGRATING THE EQUATION OF GRADUALLY VARIED FLOW” 


by Ven Te Chow 
(Proc. Paper 838) 


VEN TE cHow,! A.M. ASCE.—The writer wishes to express his deep ap- 
preciation for the constructive comments presented by the discussers. 

The major advantage of the direct method is to perform the computation of 
a flow profile through only a single operation or step. This will save a great 
deal of effort in developing a family of flow characteristic curves such as 
those demonstrated by Bakhmeteff(13)2 for canal design. The writer recog- 
nizes at the same time the advantages of the step method when applied to 
problems in which the integration method becomes evidently disadvantageous 
or fails. He therefore believes that a hydraulic engineer should learn the 
technique of both methods and know how to use the right method in a given 
problem to its best advantage. Mr. Silvester describes a step method with 
the aid of graphs, which should help the solution of a problem by the step 
method. 

For most open channels encountered in hydraulic engineering, the hydrau- 
lic exponents vary rather slowly in the design range of flow. Moreover, as 
the equation for the length of flow profile is a function of the differences be- 
tween two varied flow functions of the same hydraulic exponent, the error in- 
troduced in the result due to variation in hydraulic exponent is compensative. 
However, when the range of flow becomes large, some of the basic assump- 
tions for gradually varied flow will not be satisfactory. For example, the 
Manning coefficient of roughness changes a great deal in practically any 
channel when the variation in depth extends to a considerable range. In fact, 
the value of n for average clean sewer pipes and drain tiles, both clay and 
concrete, has been shown to increase by as much as 28 percent from depths 
equal to 1.00 to 0.25 times the diameter (D-1). For high accuracy, any 
_method assuming a constant value of n is questionable. 

When the water approaches the top of a gradually closed conduit, the hy- 
draulic exponents may change rapidly. The situation becomes further compli- 
cated due to the creation of a conjugate normal depth which will be described 
later. In such problems, the numerical integration method developed by 
Messrs. Keifer and Chu has the great advantage, because in this method the 
effect due to the variation of hydraulic exponents is included in the variation 
of the flow function. . 

One of the disadvantages of the Bakhmeteff method is the introduction of - 
an additional item o for critical slope. This item is often confused with the 
normal critical slope Sc. The connotation of this item is rather difficult, 
particularly for beginners, to conceive. Theoretically, the Bakhmeteff equa- 
tion should be identical with the author’s equation. The presence of the third 


T, Research Associate Prof. of Hydr. Eng., Graduate College, Univ. of 


Illinois, Urbana, Ill. sph ; erie 
2. Numerals in parentheses refer to corresponding items in the List o 
References in the paper and in this discussion. 
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term in the author’s equation, however, will not introduce complicated com- 
putation, but rather simple and systematic procedure as shown in the given 
example. On the other hand, the use of the item o requires additional deter- 
mination of its value even with the aid of the Kirpich graphs. 

During a visit to Europe in the summer of 1956, the writer had an oppor- 
tunity to see the original work by G. Mouret (D-2). It was a manuscript in 
mimeographed form as shown to the writer by Professor G. de Marchi, 
Director of Istituto di Idraulica e Costruzioni Idrauliche, Politecnico di 
Milano, in Milan, Italy. It is believed that Mouret’s “Hydraulique Générale” 
as mentioned by Professor Hom-ma was based on this work and published 
later in 1927 in Paris. The problem of the point of inflection is discussed on 
pages 447-458 of the manuscript. This discussion covers also another point 
of inflection which occurs theoretically on the curve for y<y¢<yn (Fig. D-1). 
This is the point which has been brought out independently for discussion in 
American literature by Gunder (D-3) in 1943. In Belgium, however, Profes- 
sor L. J. Tison of the University of Ghent has told the writer that the point of 
inflection on a theoretical backwater curve was probably first described 
analytically by Merten (D-4) in 1906, and that the possibility of a point of in- 
flection on the curve was mentioned without theoretical treatment by Boudin 
(D-5)1 at even an earlier date in 1861. The condition for an inflection point 
as given by Professor Hom-ma is obtained by placing d2x/du2 = 0. Using the 
example in the paper, where N = 3.65, M = 3.43 and yn/ye = 3.36/2.16 = 1.555, 
the writer found another point of inflection having u = 940 or y = 940yn. Ap- 
parently, this point is also beyond practical consideration. However, this 
point is on the curve for y>yn>yc, whereas the point obtained by Professor 
Hom-ma is actually for y<yc<yp: 

Owing to the peculiar hydraulic characteristics of a circular conduit or 
any conduit with a gradually closing top, a point of inflection may be shown 
to be possible at the depth of maximum discharge. Figure D-2A shows the 
variation of the normal discharge in a conduit with a gradually closing top. 

It is known that as the depth of flow increases the discharge will increase to 
the full discharge at a depth less than the full depth yp. For a circular con- 
duit, this depth yg = 0.82y, where yo is the diameter of the conduit. There- 
after, it will reach a maximum discharge at a depth ys. For a circular con- 
duit, y§ = 0.938y,. After reaching the maximum value, the discharge will 
decrease to the full discharge when the flow touches the top of the conduit. 
Within the region of y = yg and y = yo, there are two possible normal depths 
for a given discharge, namely, the lower normal depth y, and the upper or 
conjugate normal depth y,,. At maximum discharge, the two depths become 
one depth yf. It can be demonstrated analytically that four types of flow pro- 
files are possible for a given slope. Figures D-2B, D-2C, and D-2D show 
these profiles for mild and steep slopes. It should be noted that the critical 
depth in Fig. D-2D is greater than the normal depths y,;, and Yn, but its cor- 
responding lower normal depth is less than yy and y,. Consequently, the 
corresponding critical slope should be less than the normal slope, and the 
channel slope is considered as mild. 

The writer appreciates the demonstration by Mr. Harrison for extending 
the method to horizontal channels. For adverse channels, the slope of the 


1. The Sane title of this reference was supplied by Professor G. A. 
Heyndrickx of Institut Supérieur Agronomique, Départment de Hydraulique 
Agricole, Ghent, Belgium. 
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channel bottom may be taken as negative. Thus, the general differential 
equation of gradually varied flow becomes 


av 
14555) (D-1) 


SE a fe f - /yq) SY am tant 
S [ n uj+ a constan 
. leu ten" (D-2) 
in which the varied flow function can be developed as 
F(uyN)_, = S7—SB au = ph poe as 
wie S,——aH o_ 2N (D-3) 
re 1-u 1l-u 


The first integral is equal to F(u,2N) and the second integral can be easily 


shown to be (Nii) tra? oN). Therefore, 


N+1 2N (D-4) 


F(u, TD 


N)_g = F(u,2N)+ TezF (u 
ie} 


In other words, the varied flow function for adverse slopes can be computed 
directly from the varied flow functions for positive slopes. For the simplicity 
of computation, a table of the varied flow function for adverse slopes has been 


prepared as shown in Table D-1. 
The length of water surface curve between two sections 1 and 2 is 


--2 {tuy-my)- [Flapp™.,-Flo29%. i 


= (¥6/yq) aN) [F07349).g -F%p»d)_¢,} 


Example Assume a channel of adverse slope equal to 0.0016 with the same 
cross section as the example in the paper. A vertical drop controls the water 
surface at critical depth at the downstream end of the channel. Determine 
how far upstream the depth will be 10 ft. 


Solution 1. From the example yc = 2.16 ft. and y, = 3.36 ft. 
2. From 2.16 ft. to 10 ft., the average depth is 6.08 ft., and the average 
hydraulic exponents are N = 3.83 and M = 3.57. Thus, J = 3. 04. 


3. The value of (vc/yn)™ (J/N) = 0.164 
4. Computation of the distance is as follows: 
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u Vv F(u N) F(v,J) 

y N)_g. -S, 
2.16 0.643 0.573 0.621 0.549 
10.00 2.976 3.950 1.106 1.171 
Diff.: -2.333 -0.485 -0.622 


L = -2,100 [ -2.333-(-0.485)-(0.164)(-0.622) = 3,670 ft. 


Regarding the Coriolis correction factor as mentioned by Professor 
Kolupaila, the factor a can be easily added to the equation of the length of 
flow profile as follows: 


M 
L= z {(%2%9)- [F(a 9N)-F(uy N&O /y,) (3/s) [Pry )-FCv.a)}} (D-6) 
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Table D-1 The Varied Flow Function Table for Adverse Slopes 


N 
ae 2.0 2.2 204 2.6 2.8 3.0 32 34 3.6 3.8 


0.00 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 
0.02 0.020 0.020 0.020 0.020 0.020 0.020 0.020 0.020 0.020 0.020 
0.04 0.040 0.040 0.040 0.040 0.040 0.040 0.040 0.040 0.040 0.040 
0.06 0.060 0.060 0.060 0.060 0.060 0.060 0.060 0.060 0.060 0.060 
0.08 0.080 0.080 0.080 0.080 0.080 0.080 0.080 0.080 0.080 0.080 


0.10 0.099 0.100 0.100 0.100 0.100 0.100 0.100 0.100 0.100 0.100 
0.12 0.119 0.119 0.120 0.120 0.120 0.120 0.120 0.120 0.120 0.120 
0.14 04139 0.139 0.140 0.140 0.140 0.140 0.140 0.140 0.140 0.140 
0.16 0.158 0.159 0.159 0.160 0.160 0.160 0.160 0.160 0.160 0.160 
0.18 0.178 0.179 0.179 0.180 0.180 0.180 0.180 0.180 0.180 0.180 


0.20 0.197 0.198 0.199 0.199 0.200 0.200 0.200 0.200 0.200 0.200 
022 06216 0.217 0.218 0.219 04219 0.220 0.220 0.220 0.220 0.220 
O24 06234 02236 0.237 0.238 0.239 0.240 0.240 0.240 0.240 0.240 
0026 = 00253) 06255 00256 00257 00258 04259 0.259 0.260 0.260 0.260 
0.28 06272 06274 02275 02276 0.277 0.278 0.278 0.279 0.280 0.280 


0030 04291 06293 0294 02295 0.296 0.297 0.298 0.298 0.299 0.299 
0.32 0.308 0.311 0.313 0.314 0.316 0.317 0.318 0.318 0.319 0.319 
0.34 04326 04329 0.331 0.333 0.335 0.337 0.338 0.338 0.339 0.339 


0640 04380 04384 0.387 04390 04392 0.393 04394 0.395 02396 0.396 
0.42 04397 0.401 0.405 0.407 0.409 0.411 0.422 0.413 0.414 0.414 
Oe44 00414 00419 00423 00426 00429 04430 00432 06433 02434 04435 
O46 06431 06437 06440 004d 00447 00449 00451 00452 00453 00454 
O48 O6447 00453 06458 0.461 0.464 04467 02469 06471 0.472 0.473 


0.50 06463 0.470 00475 Oc479 Oe482 0485 0.487 02489 0.491 0.492 
0052 06479 0.485 O49] 06494 02499 04502 04505 06507 0.509 0.511 
0054 06494 06501 04507 06512 06516 06520 06522 06525 0.527 0.529 
056 046509 04517 02523 04528 0.533 0.537 04540 00543 0.545 0.2547 
0.58 06524 02533 0.539 06545 00550 00554 06558 06561 0.563 0.567 


0.60 06540 06548 00555 06561 00566 0.571 0.575 0.578 0.581 0.583 
Oe61L 06547 06556 06563 00569 06575 06579 0.583 0.587 0.589 0.592 
0662 06554 04563 06571 06578 0.583 0.578 0.591 0.595 0.598 0.600 
0263 04562 06571 04579 04585 02590 0.595 0.599 0.603 0.607 0.609 
0.64 04569 0.579 0.586 0.592 0.598 0.602 0.607 0.611 0.615 0.618 


0.65 04576 04585 0.592 0.599 0.606 0.610 0.615 0.619 0.623 0.626 
0.66 0.583 0.593 0.600 0.607 0.613 0.618 0.622 0.626 0.630 0.634 
0.67 0.590 04599 0.607 0.614 06621 0.626 0.631 0.635 0.639 0.643 
0.68 0.597 0.607 0.615 0.622 0.628 0.634 0.639 0.643 0.647 0.651 
0.69 0.603 0.613 0.621 0.629 0.635 0.641 0.646 0.651 0.655 0.659 
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Table D-2 The Varied Flow Function Table for Adverse Slopes 


1-000 0.785 0.797 0.808 0.818 0.826 0.834 0.842 0.849 0.856 0.862 
1.005 0.788 0.799 0.810 0.820 0.829 0.837 0.845 0.852 0.858 0.864 
1.010 0.790 0.801 0.812 0.822 0.831 0.840 0.847 0.855 0.861 0.867 
1.015 0.793 0.804 0.815 0.824 0.833 0.843 0.850 0.858 0.864 0.870 
1.020 0.795 0.807 0.818 0.88 0.837 0.845 0.853 0.860 0.866 0.872 


1.03 0.800 0.811 0.822 0.832 0.841 0.850 0.857 0.864 0.871 0.877 
1.04 0.805 0.816 0.829 0.837 0.846 0.855 0.862 0.870 0.877 0.883 
1.05 0.810 0.821 0.831 0.841 0.851 0.859 0.867 0.874 0.881 0.887 
1.06 0.815 0.826 0.837 0846 0.855 0.864 0.871 0.879 0.885 0.891 
1.07 0.819 0.831 0.841 0.851 0.860 0.869 0.876 0.883 0.889 0.896 


1.08 0.824 0.836 0.846 0.856 0.865 0.873 0.880 0.887 0.893 0.900 
1.09 0,828 0.840 0.851 0.860 0.870 0.877 0.885 0.892 0.898 0.904 
1.10 0.833 0.845 0.855 0.865 0.874 0.881 0.890 0.897 0.903 0.908 
1-11 0.837 0.849 0.860 0.870 0.878 0.886 0.894 0.900 0.907 0.912 
Lee 0.842 0.854 0.864 0.873 0.882 0.891 0.897 0.904 0.910 0.916 


1.13 0.846 0.858 0.868 0.873 0.886 0.895 0.902 0.908 0.914 0.919 
114 0.851 0.861 0.872 0.881 0.890 0.899 0.905 0.912 0.918 0.923 
1015 0.855 0.866 0.876 0.886 0.895 0.903 0.910 0.916 0.922 0.928 
1.16 0.859 0.870 0.880 0.890 0.899 0.907 0.914 0.920 0.926 0.931 
1.17 0.864 0.874 0.884 0.893 0.902 0.911 0.917 0.923 0.930 0.934 


1.18 0.868 0.878 0.888 0.897 0.906 0.915 0.921 0.927 0.933 0.939 
1.20 0.876 0.886 0.896 0.904 0.913 0.921 0.928 0.934 0.940 0.945 
1.22 0.880 0.891 0.900 0.909 0.917 0.929 0.932 0.938 0.944 0.949 
1.24 0.888 0.898 0.908 0.917 0.925 0.935 0.940 0.945 0.950 0.955 


1.26 0.900 0.910 0.919 0.927 0.935 0.942 0.948 0.954 0.960 0.964 
1.28 0.908 0.917 0.926 0.934 0.945 0.948 0.954 0.960 0.965 0.970 
1.30 0.915 0.925 0.933 0.941 0.948 0.955 0.961 0.966 0.981 0.975 
1632 0.922 0.931 0.940 0.948 0.955 0.961 0.967 0.972 0.976 0.980 
1.34 0.930 0.939 0.948 0.955 0.962 0.967 0.973 0.978 0.982 0.986 


1.36 0.937 0.946 0.954 0.961 0.968 0.973 0.979 0.983 0.987 0.991 
1.38 0.944 0.952 0.960 0.967 0.974 0.979 0.985 0.989 0.993 0.996 
1.40 06951 0.959 0.966 0.973 0.979 0.984 0.989 0.993 0.997" 1.000 
1.42 04957 0.965 0.972 0.979 0.984 0.989 0.995 0.998 1.001 1.004 
144 0.964 0.972 0.979 0.984 0.990 0.995 1.000 1.003 1.006 1.009 
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Table D-2 The Varied Flow Function Table for Adverse Slopes 


N 
mo EO” £2 --GL5, G0 58 


0.00 0.000 0.000 0.000 0.000 0.000 
0.02 0.020 0.020 0.020 0.020 0.020 
0.04 0.040 0.040 0.040 0.040 0.040 
0.06 0.060 0.060 0.060 0.060 0.060 
0.08 0.080 0.080 0.080 0.080 0.080 


0.10 0.100 0.100 0.100 0.100 0.100 
0.12 0.220 0.120 0.120 0.120 0.120 
0.14 0.140 0.140 0.140 0.140 0.140 
0.16 0.160 0.160 0.160 0.160 0.160 
0.18 0.180 0.180 0.180 0.180 0.180 


0.20 06200 0.200 0.200 0.200 0.200 
0622 0.220 04220 04220 0.220 0.220 
024 04240 0.240 0.240 0.240 0.240 
0.26 0.260 0.260 0.260 0.260 0.260 
0.28 0.280 0.280 0.280 0.280 0.280 


0.30 04300 0.300 0.300 0.300 0.300 
0.32 046320 04320 04320 0.320 0.320 
0.34 02339 0.340 04340 0-340 0.340 
0.36 0.359 0.360 0.360 0.360 0.360 
0.38 0.378 0.379 0.380 0.380 0.380 


0.40 04397 0.398 0.398 0.400 0.400 
0.42 0.417 0.418 0.418 0.419 0.2420 
0.44 02436 06437 02437 0.439 0.440 
0046 046455 02456 06457 04458 04459 
048 06474 06475 06476 0.478 0.4479 


0.50 0.493 0494 02495 0-497 0.498 
0.52 0.512 0.513 0.515 0.517 0.518 
0054 046531 04532 02533 04536 0.537 
0.58 0.567 04569 0.570 0.574 0.576 


0.60 0.585 0.587 0.589 0.593 0.595 
0.62 0.603 0.605 0.607 0.611 0.613 
0.63 0.612 0.615 0.616 0.620 0.622 
0.64 0.620 0.623 0.625 0.629 0.631 


0.65 0.629 0.632 0.634 0.638 0.640 
0.66 0.637 0.640 0.643 0.647 0.750 
0-67 06646 06649 0.652 0.656 0.659 
0.68 0.654 0.657 0.660 0.665 0.668 
069 0.662 0.665 0.668 0.674 0.677 


ASCE 


DISCUSSION 


Table D-2 The Varied Flow Function Table for Adverse Slopes 


N 
ANS 4.0 Pro ies Lee 5.0 5.5 


0.70 0.670 0.673 04677 0.682 0.686 
0.71 0.678 0.681 0.685 0.690 0.694 
0.72 0.686 0.689 0.694 0.699 0.703 
0.73 02.694 0.698 0.702 0.707 0.7122 
0.74 0-702 0.705 0.710 0.716 0.720 


0.75 0.709 0.712 0.717 0.724 0.728 
0.76 0.717 0.720 0.725 0.731 0.736 
enim eOereL- Onrer 02723 02739 0-744 
0.78 0.731 06735 04740 06747 04752 
079 06738 06742 0.748 06754 0.760 


0.80 04746 02750 0.755 0.762 0.768 
0.81 0.4753 06757 0.762 04770 0.776 
0.82 0.760 0.764 0.769 0.777 0-783 
0.83 0.766 04771 0.776 0.784 0-79 
0.84 02773 00778 0.783 0.791 0.798 


0.85 0.780 0.784 0.790 0.798 0.805 
0.86 0.786 0791 0.797 0.804 0.812 
0.87 0.793 0.797 0.803 0.811 0 2819 
0.88 0.799 0.803 0.810 0.818 0 2826 
0.89 0.805 0.810 0.816 0.825 0 2832 


(@) e7 0 e811 0 2816 10) 2822 0 2831 10) 2839 
0.91 0.817 0.821 0.822 0.837 0 0845 
0.92 0.823 0.828 0.834 0.844 0.851 
0.93 04829 0.833 0.840 0 0850 02857 
0.94 02835 04.840 0 0846 0.856 0-864 


0.950 0.840 0.845 02852 0.861 0.869 
0.960 0.846 0.861 0.857 0.867 0.875 
0.970 0.851 0.866 0.863 0.972 0.881 
0.975 0854 04859 0.866 0.875 0.883 
0.980 02857 0.861 0.868 0.878 0.886 


0.985 0859 0.863 0.870 0.880 0.889 
0.990 0.861 0.867 0.873 0.883 0.891 
0.995 04864 0.869 0.876 0.885 0-894 
1.000 0.867 0.873 0.879 0.887 0.897 
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Table D-2 The Varied Flow Function Table for Adverse Slopes 


N 
a LO a ate) ha 50 5.5 


1.000 0.867 0.873 0.879 0.887 0.897 
1.005 0.870 0.874 0.881 0.890 0.899 
1.010 0.873 0.878 0.884 0.893 0.902 
1.015 0.875 0.880 0.886 0.896 0.904 
1.020 0.877 0.883 0.889 0.898 0.907 


1.03 0.88 0.887 0.893 0.902 0.911 
1.04 0.888 0.893 0.898 0.907 0.916 
1.05 0.892 0.897 0.903 0.911 0.920 
1.06 0.896 0.901 0.907 0.915 0.924 
1.07 0.901 0.906 0.911 0.919 0.928 


1.08 0.905 0.910 0.916 0.923 0.932 
1.09 0.909 0.914 0.920 0.927 0.936 
1.10 0.913 0.918 0.923 0.931 0.940 
1.11 0.917 0.921 0.927 0.935 0.944 
1.12 0.921 0.926 0.931 0.939 0.948 


1.13 0.925 04929 0.935 0.943 0.951 
1.14 0.928 0.933 0.938 0.947 0.954 
1.15 0.932 0.936 0.942 0.950 0.957 
1.16 0.936 0.941 0.945 0.953 0.960 
1.17 0-939 0.944 0.948 0.957 0.963 


1.18 0.943 0.947 0.951 0.960 0.965 
1.19 0.947 0.950 0.954 0.963 0.968 
1.22 0.956 0.957 0.964 0.972 0.976 
1.24 0.962 0.962 0.970 0.977 0.981 


1-26 0.968 0.971 0.975 0.982 0.986 
1.28 0.974 0.977 0.981 0.987 0.990 
1.30 0.979 0.978 0.985 0.991 0.994 
1.22 0.985 0.986 0.990 0.995 0.997 
1.34 0.990 0.992 0.995 0.999 1.001 


1.36 0.994 0.996 0.999 1.002 1.005 
1.38 0.998 1.000 1.003 1.006 1.008 
1.40 1.001 1.004 1.006 1.009 1.011 
1-42 1.005, 1.008 2.010 1.012 1.014 
1-44 1.009 1.013 1.014 1.016 1.016 
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FIG.D-2 FLOW PROFILES OF A CLOSED CONDUIT. 
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Discussion of 
“FREE OUTLET AND SELF-PRIMING ACTION OF CULVERTS” 


by Wen-Hsiung Li and Calvin C. Patterson 
(Proc. Paper 1009) 


_ J. E. FLACK,* J.M. ASCE.—The authors’ presentation deals with some 
jaspects of culvert flow that are not ordinarily considered when designing 
these structures.(1) Experiments with which the writer is familiar, conduct- 
ed at the Iowa Institute of Hydraulic Research by graduate students in the 
Department of Mechanics and Hydraulics, are hereafter described because of 
the similarity and extension to the results as presented in the paper under 
discussion. 


Piezometric Grade Line 


The Iowa experiments determined the location of the intersection of the 
hydraulic grade line with the plane of the culvert outlet, for free flow, for a 
4 inch round conduit by Rueda(2) and for a 3 inch square box culvert model by 
the writer.(3) The evaluation of the factor m, i.e., the ratio, expressed as a 
decimal, of the height of the intersection of the hydraulic grade line with the 
plane of the outlet to the total height of the pipe or culvert, was determined 
by extending the hydraulic grade line to the plane of the outlet. The hydraulic 
grade line in the culvert barrel was determined by piezometer connections at 
points along the length of the culvert. All measurements were made with the 
culvert horizontal and flowing full. 

The variation of m with the Froude number of the flow was determined for 
both the round conduit and the box culvert, and the results are shown on Fig- 
ure 1. Also shown are two values for m as determined by the staff of the 
Institute using the relaxation method for a two-dimensional, frictionless, con- 
duit outlet. The authors’ curves are also reproduced for comparison 
purposes. 

_ For the round pipe, the authors’ curve is lower than that of Rueda at 
Froude numbers greater than 1.0. For the box culvert, the departure of the 
authors’ curve from the plotted values as obtained in the Iowa experiments 
is significant above Froude numbers of 1.5, with the Iowa results consistently 
lower. 

Since the experiments were similar the discrepancy in results can, per- 
1aps, be attributed to the different methods of evaluating m, although some of 
he authors’ plotted values are quite close to those as obtained in the Iowa 
‘xxperiments. Further experimentation is suggested to determine values of m 
or higher values of the Froude number of the flow before the authors’ curves 
ire used to “furnish proper values of m for design.” 


« Instr., Dept. of Civ. Eng., Univ. of Colorado, Boulder, Colo. 
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Fig. 2 Vortex Over a Well-Rounded Inlet. 
(Culvert Flowing Full) 
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Priming 


Self-priming of culverts was also noted in the Iowa experiments. A sig- 
nificant phenomena, apparently not encountered by the authors, was the 
formation of vortices over the culvert inlet. These had the effect of counter 
acting the priming action by feeding air to the zone of separation. In most 
cases, the Iowa models differed from those of the authors only in that the 
slope of the top of the wing-walls was continued above the culvert much like 
a highway embankment. Figure 2 is a photograph of a vortex over a well 
founded inlet of a culvert flowing full. 

A model of a standard culvert inlet, with square entrance and flaring 
wingwalls, was modified and the results can be summarized as follows: 


1. The standard inlet will perform more efficiently with the addition of a 
smooth well-rounded entrance. It is more important to round the top than 
the sides. 

2. A well-rounded inlet, with rounding on the top and sides, and a headwal 
as high as the design head and approximately twice as wide as the culvert wil 
give the most efficient performance, and vortex action will be eliminated. 


Another aspect noted was the scale effect in evaluating prototype per- 
formance for culverts. In the Iowa tests a 6 inch model box culvert, as 
tested by Maksoud,(4) gave appreciably better performance than a 3 inch 
model of similar design. It might be expected that results of similar signifi- 
cance can be expected with larger models and prototype structures. 
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JOHN L. FRENCH,* A.M. ASCE.—The authors are to be congratulated on 
their contributions to culvert hydraulics. These contributions are two-fold: 
first, they have devised an ingenious experimental method for determining tl 
effective pressure at the culvert outlet and secondly, they have contributed 
significantly to a better understanding of the various complex phenomena by 
which the regime of flow in a culvert may change from part full to full condu 
flow. 

During the past two years the Bureau of Public Roads has sponsored an i 
vestigation of pipe culvert flow in the hydraulic laboratory of the National 
Bureau of Standards. While the project is not yet completed, some of the 


experimental data now available are pertinent to the problems investigated | 
by the authors. 


* Hydr. Engr., National Bureau of Standards, U. S. Dept. of Commerce, 
Washington, D. C. 
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In a sufficiently long culvert, flowing full, three distinct zones of flow 
Jong the culvert may be observed. The first may be termed the zone of flow 
stablishment in which boundary layer development plays an important role, 
7ith consequent variation of the velocity profile with distance from the inlet. 
1S the flow proceeds downstream the developing boundary layer grows pro- 
ressively thicker and eventually reaches the center of the pipe. At this 
oint the velocity profile becomes stable with no further variation of shape 
vith distance from the inlet. This constancy of the velocity profile is a 
haracteristic of the second zone of flow, sometimes termed the zone of 
istablished flow. 

At the outlet, assuming a non-submerged condition with either free fall or 
nm apron support, non-hydrostatic pressure distribution prevails. Owing to 
he necessary redistribution of pressures required by the change from hydro- 
tatic pressure distribution in the zone of established flow to the non- 
ydrostatic pressure distribution at the outlet, a terminal length exists near 
he outlet in which accelerative effects become important with consequent 
istortion of the axisymmetrical velocity profile of the zone of established 
low to a non-symmetrical velocity distribution at the outlet. It is the phe- 
omenon occurring in the terminal length with regard to redistribution of 
ressures and velocities which makes necessary the determination of an 
ffective pressure at the outlet for use with the one-dimensional method of 
nalysis of full conduit flow in culverts. The authors have investigated this 
roblem in an ingenious and simple experimental manner for a conduit 20 
iameters in length. 

The authors’ basic equation for defining m, is 


7A 
V (15) 
Hy ok = On Qe 


+ mD + Ay 

hich follows directly from the authors’ equation (1) on the assumption that 
1e approach velocity is negligible. In equation (15) the energy correction 
ctor @, applies to the outlet section. Under these circumstances the 
uthors correctly define mD as the mean piezometric head at the outlet sec- 
on; i.e. 


1 P41 
mD = 6 \ = +z) dQ, (16) 
Q 


here p, is the pressure at the outlet section a vertical distance z above the 
wert. At the same time the authors state the “. . . effective value of mD for 
mg culverts can also be obtained by projecting the hydraulic grade line to 
1e outlet section.” There would appear to be a conflict of ideas here. 

Let BD be the elevation above the invert at the outlet section of the extend- 
1 straight line hydraulic gradient of the zone of established flow. It then 
lows that 


{, vy? (17) 
2g 
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where ({ is the length of the terminal segment and h is the piezometric head 
at the upstream end of the terminal length, the datum being taken as the 
invert of the outlet section. 

For the segment of the culvert including only the zone of flow establish- 
ment and the zone of established flow it is clear that 


2 Urakg) sing? (18) 


H +SL-h=h, +a", ek ee Wa Qe? 


1 


where a eT is the kinetic energy factor for the zone of established flow and he 
is the so-called entrance loss. Eliminating h through equation (17), 


2 2 (19) 
V Ean) 
H, + SL =a' 1 2g +BD+ th, +£5 90) 


Subtracting equation (15), 


2 2 
ks 1 ; EY (20) 
Da BD es (a, - a3) - Be = (hot fe | 


D 2g 


and it is clear that m and # are not identical. The distinction between the 
two is that mD is the mean piezometric head at the outlet, as defined by 
equation (16) which follows from equation (15), while 8 D is an effective 
piezometric head defined by 


p 2 ' 2 
pD=% | Gisnags¥ (ay -o' + | 4-0, +f BE | (21) 
Q 


which follows from equations (16) and (20) through the definition of the effec- 
tive piezometric head at the outlet section as being the elevation above the it 
vert of the point at which the extended straight line piezometric gradient of 
the zone of established flow pierces the outlet plane. 

Consideration of equation (21) will show that the experimental problem of 
determining BD, as well as mD, of course, stems entirely from the deviatiot 
of the pressure distribution from hydrostatic at the outlet. That is, if hydro 
static pressure distribution prevails throughout the terminal length, vertical 
accelerations in the terminal length will be absent and it follows that curva- 
ture of the stream lines in the terminal length can not occur. Hence the axi 
symmetrical velocity distribution of the zone of established flow will continu 
undistorted through the terminal length, and referring to equation (21) 


1 Z 1 
Q | $a0-5, 
Q 
’ 
Oy -a, = 0, 
and Lv? 
Ay (h, ars 0, 
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from which 


p 
ial 1 


p 
If it is further assumed that LU | ahs. dQ is zero, B = 1/2, the value 
Q Q pgD 
)}commonly used in the conventional one-dimensional pipe flow theory for a 
non-submerged outlet. 

If it is assumed that the outlet is submerged by a tail water depth T and 
that hydrostatic pressure distribution occurs at the outlet, it follows from 
equation (22) that B= T/D the value commonly assumed for a submerged 
outlet. 

With reference to the experimental method used by the authors let equa- 
tion (19) as it refers to a particular type of outlet condition be 


2 
AAW. LV (23) 
H,+SL=o,55+ Petes hon) 
and with a second type of outlet condition 
2 2 
eV: LV 
Ho i: Slaalaige ReaD Meetings) s 


Subtracting equation (23) and rearranging 


2 
V 1 ’ 
B,D = By D - (Hy - Hy) +55 (9 -O 4) + (yy - hey). 


If the culvert is sufficiently long to have a zone of established flow, o& 9 = 
ay and 


B,D = B,D - (Hy - Hy) + (hg, - hy). (24) 


If it is now assumed that the subscript 2 refers to a submerged outlet in 
which the tail water elevation above the invert is T, equation (24) becomes on 
the assumptions that (a) 8 gD = T and (b) heg = hey, 


B,D aed We (Hy iS H,), 


the relation used by the authors to compute mD for the various outlet 
conditions. ine 
Since both of these assumptions were made by the authors in the derivation 
of equation (4) as well as the assumption that the velocity distribution at a 
outlet face is identical for the submerged and nonsubmerged cases, it wou 
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appear that the authors novel experimental method would determine f with 
perhaps greater accuracy than m. 

The authors have used a conduit length of 20 diameters in determining the 
effective pressure at the outlet face. In view of the work of Morris” it would 
appear questionable if the fully developed turbulent velocity profile charac- 
teristic of the zone of established flow actually existed at the upstream end of 
the terminal length. Some influence of velocity distribution at the upstream 
end of the terminal length on B or m would be expected. Whether the magni- 
tude of this effect would be such as to have practical significance is not clear. 

In this connection it may be observed that in the above discussion of f it 
was assumed that the culvert barrel was sufficiently long so that a zone of 
established flow existed. It is obvious that for culverts so short as not to 
possess a zone of established flow that the concept of the effective pressure 
head at the outlet as being the elevation above the invert at which the extended 
straight line piezometric gradient of the zone of established flow pierces the 
plane of the outlet is without direct physical significance. For such short cul- 
verts equation (19) employs a resistance coefficient f, which is no longer 
directly applicable to any segment of the culvert barrel, and it employs a 
kinetic energy factor oy and an entrance head loss he, neither of which are 
now invariant with culvert length. Under these circumstances the utility of 
equation (19) with reference to short barrels is open to question, as is also, 
of course, the physical significance of equation (17) which defines B in terms 
of the resistance coefficient of the zone of established flow. 

The authors describe three types of filling or priming phenomena charac- 
terizing the change from part full sluice type flow to full conduit flow. In 
tests made by the writer with a smooth lucite model culvert 5.5 inches in 
diameter and 63.9 diameters long, two additional types of filling phenomena 
were observed. These consisted of: 


A. Filling through direct backwater effect. The authors under their type 1 
priming phenomena for mild slopes describe the filling of a model culvert by 
means of a hydraulic jump. The presence of a jump in the authors’ tests with 
sub-critical slopes would imply supercritical flow in the upstream portion of 
the barrel, with the presence of a typical vena contracta near the inlet. This 
type of flow may for convenience be termed contracted flow at the inlet. In 
the writer’s tests for mild slopes this filling process by means of a hydraulic 
jump was also observed but as the slope was progressively decreased a point 
was reached where contracted flow at the inlet and hence a hydraulic jump 
did not occur and the filling mechanism consisted simply of a backwater ef- 
fect in which the water surface near the inlet approached the barrel crown as 


@/p! a was increased. 

It was observed experimentally that the filling mechanism consisted of the 
non-contracted flow backwater effect described above at zero slope. Ata 
slope of 0.25% contracted flow occurred with filling by means of a hydraulic 
jump. The criterion which determines if filling by this method will occur is 
apparently the depth of water, due to backwater effect, at the section where 
the vena contracta would ordinarily form, when the inlet becomes submerged. 
If this depth due to backwater is greater than the conjugate or alternate depth 


7, Preliminary Flow Tests on a Model Culvert, by Henry M. Morris, Project 
Report No. 7, St. Anthony Falls Hydraulic Laboratory, Univ. of Minn. 
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corresponding to the rate of discharge and the vena contracta depth, a jump 
cannot form because of momentum considerations, hence contracted flow 
cannot occur, and the filling mechanism consists simply of the backwater 
curve approaching the culvert crown at the inlet as the rate of discharge is 
)sufficiently increased. However, if the depth referred to is less than the 
iconjugate depth then a jump occurs, with contracted flow at the inlet, for suf- 
ficiently long culverts on sub-critical slopes, and filling occurs when the rate 
of discharge is increased sufficiently to cause the water surface immediately 
downstream of the jump to approach the culvert crown. 


B. Combined free surface and full conduit flow. As the barrel slope is in- 
creased through critical slope, sub-critical flow can no longer exist in the 
downstream portion of the barrel, and hence filling by means of a hydraulic 
jump can no longer occur. For slopes in the range 0.6% to 1.25% the filling 
mechanism, for the 63.9 diameter long conduit used by the writer, was ob- 
served to consist of the following: The flow was supercritical throughout the 
culvert barrel, and because of boundary roughness the maximum depth of 
flow occurs at or slightly upstream from the outlet. Also owing to the super- 
critical flow, the water surface is subject to standing shock waves generated 
by the channel contraction at the inlet. As the flow is increased a point is 
reached where initial contact of the water surface with the culvert crown 
occurs one or two diameters upstream from the outlet. The contact point 
advances a short distance upstream and becomes stationary with the portion 
of the culvert barrel downstream of the contact point flowing full. This type 
of flow was characterized by entrained air being passed through the full con- 
duit flow and the presence of an intermittent vortex in the upstream pool 
feeding air into the culvert barrel and thence to the air pocket extending from 
the inlet to the contact point marking the upstream end of the segment of the 
barrel flowing full. The pressure in the air pocket was appreciably sub- 
atmospheric. Flow conditions thus consisted of the portion of the barrel im- 
mediately above the outlet operating in full conduit flow and an upstream 
portion of the barrel operating in free surface flow, with the air space above 
the free surface under a non-steady, sub-atmospheric pressure with air be- 
ing withdrawn at the downstream end of the 50 or 60 diameter long air pocket 
by air entrainment through the full conduit portion of the barrel and an inter- 
mittent supply of air coming to the air pocket from vortex action in the up- 
stream pool. 

After initial filling at the outlet, as the rate of discharge was increased the 
contact point of the free water surface with the culvert crown moved pro- 
gressively upstream until it reached a distance of 15 or 20 diameters from 
the outlet. Any further increase in the rate of discharge then caused the cul- 
vert to fill throughout its entire length. : 

The filling mechanism in this case is connected with the internal forces 
controlling the relative magnitude of the rate of discharge Q, in the free 
surface segment of the barrel, and the rate of discharge Q2; of the full con- 
juit flow segment. If for a given location of the contact point of the free sur- 
lace with the culvert crown (marking the upstream end of the full conduit flow 
segment), Q; > Q2 the contact point will move upstream; if Q; a Q2 ay CAs 
ibrium condition will have been established and the contact point will be sta- 
ionary, and if Qo > Qy the contact point will move downstream. The dis- 
charge of the full conduit segment will of course be a function of the net 
vharge Qo & 5 

i full conduit flow segment. That 
ressures acting to force water through the 
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is, the effective pressure on the upstream end of the full conduit section and 
hence Q» will be determined by the air pressure over the free surface seg- 
ment and by momentum considerations relating to the depth and velocity of 
the free surface flow immediately upstream from the full conduit flow seg- 
ment. In like manner the rate of discharge Q, of the free surface flow will 
be determined, for a given inlet, by the depth of submergence of the inlet and 
the air pressure over the free surface flow. 

With regard to the effect of the sub-atmospheric pressure over the free 
surface, consider the typical case with this type of filling mechanism in which 
a short segment of the barrel immediately upstream from the outlet is flow- 
ing full and the rest of the barrel is flowing with a free surface. Assume 
first that the air pressure over the free surface flow is atmospheric and that 
the flow is in equilibrium. That is, Qy = Qg and the contact point of the free 
surface with the barrel crown is stationary. If it is now assumed that a sub- 
atmospheric pressure develops over the free surface, it is obvious that Q, 
will increase and Q2 will decrease. That is a non-steady flow condition will 
have developed in which Qy > Qy and consequently the contact point will move 
upstream either to a new equilibrium point, the internal forces permitting, or 
regime change to full conduit flow over the whole length of the barrel will 
occur. It is apparent that a decrease in pressure over the free surface will 
have the effect of causing the contact point to move upstream and hence would 
be expected to promote change to full conduit flow at lower inlet submer- 
gences. 

That the air pressure over the free surface appreciably affects the filling 
mechanism was demonstrated experimentally by installing a small air vent 
over the free water surface, With the vent open for the passage of air from 
the atmosphere to the air pocket, filling occurred for a slope of 1.25% at an 
inlet submergence, 0.65D higher than when the vent was closed. Measure- 
ment of these sub-atmospheric pressures over the free surface indicate that 
they are not of inconsequential magnitude; the pressure head, referred to at- 
mospheric, varying from -0.18D for a slope of 0.6% to -0.57D for a 1.25% 
slope. These measurements were made immediately after initial contact of 
the water surface with the barrel crown had occurred near the outlet, and 
represent an equilibrium condition in which the contact point was stationary 
(or more precisely, oscillating about a stationary point) approximately four 
to six diameters upstream from the outlet. They do not represent the pres- 
sures just prior to regime change to full conduit flow. That such pressures 
would show a further decrease is indicated by a preliminary set of data for 
the 0.6% slope in which the pressure head decreased from -0.18D for initial 
contact of the water surface with the crown to -0.25D just prior to regime 
change to full conduit flow throughout the culvert. 

These sub-atmospheric pressures over the free surface flow are, as has 
been noted, the result of pumping air out of the air pocket through air en- 
trainment at the upstream end of the full conduit flow section and an inflow 
of air through vortex action at the inlet. 

It is to be noted in the above description that three conditions are neces- 
sary for this particular type of filling mechanism to occur. First the barrel 
length, slope and roughness must be such that initial contact of the water 
surface with the crown occurs at or slightly upstream from the outlet sec- 
tion. Second, the rate of discharge required to cause initial filling at the out- 
let section must be such that the internal forces corresponding to this rate of 
discharge may reach an equilibrium point where Q; = Qg for a contact point 
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less than L distant from the outlet. And third, as the discharge is increased 
a point must be reached where Qy > Qg throughout the culvert, with conse- 

, quent regime change to full conduit flow. 

| For the 63.9 diameter long, smooth lucite culvert used in the writer’s 

_ tests, regime change to full conduit flow by the phenomenon described in B 
above failed to occur as the slope was increased above 1.25% through inability 
of the flow to satisfy the second of the three requirements for type B filling. 
That is, regime change to full conduit flow throughout the entire barrel 
occurred immediately after initial filling of the conduit at the outlet section, 
in the same manner as described by the authors under “Self-priming with a 
divergent flow.” In the tests made by the writer the initial section of filling 
was invariably within 2 to 4 diameters of the outlet. 

In the writer’s tests piezometer openings were placed on the culvert crown 
at points 0.7D, 7.6D, 24.2D and 43.2D from the inlet. These piezometer open- 
ings made it possible to measure the air pressure within the conduit during 
part full flow, as well as during the filling operation described under B 
above. Preliminary analysis of the data indicates that the air pressure, for 
a square edge inlet, may be represented to an approximation by a relation of 
the form 


: 


oe 


= Q 
gD fw ale Sp gh 


a) 


where @ is for the single 63.9 diameter long smooth culvert used, a function 
only of the distance from the inlet, and pg is the air pressure over the free 
surface, referred to atmospheric as a datum. 

The experimentally determined values of 0, at distances of 0.7D and 7.6D 
from the face of the inlet were 3.8 x 10-7 and 0.96 x 10-7 respectively, indi- 
cating first that a large adverse air pressure gradient exists near the inlet 
and secondly that the sub-atmospheric pressure head near the inlet is sub- 
stantial for the higher rates of flow; amounting to -0.25D at the first piezome- 


ter opening for Q/p*/ oe 9.5, a value of q/p*/ 2 still within the range of the 
outlet control filling mechanism described in the authors’ type 2 filling. 

In the writer’s tests with a longer and greater diameter barrel than used 

‘by the authors the precise filling mechanism described by the authors under 
“(3) Self-priming due to standing surface waves,” in which priming occurred 
near the inlet due to a wave crest approaching the barrel crown, was not ob- 
served. In this connection it may be noted that in the writer’s tests the 
water surface for a distance of 6 or 8 diameters downstream from the inlet, 
except for a segment immediately downstream from the inlet of approximate 
length 1.5D, could not be observed visually because of a sheet of water flow- 
ing around the pipe boundary. 

Since the rate of discharge required to cause the water surface to ap- 
proach the outlet crown increases as the slope is increased, it is apparent, 
as was observed experimentally, that the inlet section must be subjected to 
increasingly large sub-atmospheric pressures as the slope increases with 
type 2 filling. As observed in the writer’s tests with a longer barrel than 
the authors’, this rapid increase in sub-atmospheric pressure in the inlet 


‘section as @/p*/ z approaches and increases above 9.5 appears to be an im- 
portant part of the internal mechanism operating within the culvert to cause, 
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as the slope is increased above 1.8%,a change in the filling mechanism from 
outlet control in which filling occurs through the water surface approaching 
the culvert crown to a complex inlet filling mechanism involving both 
pneumatic and hydraulic phenomena. 

For the 63.9D long smooth culvert used in the writer’s tests, initial filling 
at the inlet occurred at slopes of 2% and greater. Owing to the supercritical 
flow, the water surface was subjected to shock waves caused by the channel 
contraction at the inlet, as observed by the authors. These shock waves ex- 
hibited the usual pattern of such waves, with the water surface transverse to 
the barrel axis deviating appreciably from a horizontal line. The waves were 
of appreciable magnitude near the inlet and decreased in amplitude as the 
outlet was approached. The barrel segment in the region from the inlet toa 
section approximately 8 diameters downstream was observed as previously 
stated to be under sub-atmospheric air pressures of relatively large magni- 
tude just prior to regime change from part full to full conduit flow. Violent 
entrainment of air was observed near the inlet as reported by the authors. At 
a section 8 diameters downstream from the inlet the water was observed to 
contain many air bubbles. This relatively large volume of entrained air ob- 
served at this section was for the most part discharged from the water before 
a section approximately 12 to 15 diameters downstream from the inlet was 
reached. 

Inasmuch as an air-carrying vortex is in general not present at the high 
pool levels required to cause regime change on the slopes now considered, 
the presence of entrained air in the inlet region indicates a flow of air above 
the water surface in the upstream direction. There would thus appear to be a 
relatively large local circulation of air in the barrel segment extending from 
the inlet to a section 12 to 15 diameters downstream. That is, air is en- 
trained by the water near the inlet and carried 12 to 15 diameters downstrear 
where it enters the air space above the free surface and then flows upstream 
to the inlet region where it is re-entrained. 

As the rate of discharge is increased a point is eventually reached where 
the sluice type flow abruptly changes to full conduit flow. For all slopes 
tested of 2% to 9.39%, the initial change to full conduit flow occurred 6 or 8 
diameters downstream from the inlet. Invariably after the initiation of full 
conduit flow at this section the region of full conduit flow spread abruptly up- 
stream and downstream filling the entire conduit. Also invariably, the onset 
of full conduit flow was characterized by the abrupt appearance of foamy 
water throughout the conduit cross section. As indicated previously in no 

_ case was the filling mechanism for the 63.9 long barrel used by the writer 
observed to be the result of simple wave action as described by the authors 
for their shorter barrels under “Self-priming due to standing surface waves” 


with filling being caused by the gradual approach, as @/p*/ : was increased, 
of a wave crest to the culvert crown. 

That the relatively large adverse air pressure gradient in the 5 to 10 
diameter length segment immediately downstream from the inlet exerts a 
marked influence on the filling mechanism for this type of flow was demon- 
strated by allowing air to flow into the region from the atmosphere through ~ 
the piezometer openings on the top of the culvert. This effect was first ; 
demonstrated by Straub, Anderson and Bowers,® who showed that for a Squar 


8. Effect of Inlet Capacity of Culverts on Steep Slopes by L. G. Straub, 
(footnote contd. pg. 3! 
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€dge inlet to a smooth lucite circular barrel of 105 diameters in length on a 


; 

| 4% slope, full conduit flow resulted at @/p°/? = 9.35 for normal operation, 

| while it did not when the culvert was “ventilated,” presumably by allowing 

air to flow into the sub-atmospheric pressure region above the free surface 

| through piezometer openings in the crown of the barrel. 

This demonstrated effect of the adverse air pressure gradient near the in- 

_ let, in lowering the pool level at which filling occurs, would be qualitatively 

expected from considerations of the effect of such a gradient on the down- 
stream flow. In this connection it is obvious from energy considerations 
that for a given rate of flow, the effect of the adverse pressure gradient over 
the free surface will be to decrease the specific energy of the flow, and hence 
increase the depth of flow at points in the culvert barrel downstream of the 
inlet as compared with the depths that would have prevailed at the same rate 
of flow, had the adverse pressure gradient not been present. The increase in 
depth, for a given discharge, caused by the adverse pressure gradient at a 
given downstream point would be expected, of course, to cause regime change 


to full conduit flow at a lower value of Q/ p/ (and H/D) than would be the 
case had the gradient not been present. 

As observed in the writer’s tests, the precise mechanism by which filling 
occurs near the inlet, involving as it does an adverse air pressure gradient 
of appreciable magnitude, a relative large local air circulation in the barrel 
for some 12 to 15 diameters downstream, and an air entrainment process 
which involves the energy losses associated with the entrainment and trans- 
port of air from a region of low pressure to a region of higher pressure, well 
deserves further experimental study. 

Filling of this type might be expected to be relatively insensitive to length 
except for barrel lengths which placed the outlet within the 12 to 15 diameter 
long segment of conduit being subject to the observed severe local circulation 
of air above the free water surface. Reduction of barrel length sufficiently 
to place the outlet in this region would be expected to lower the adverse 
pressure gradient and hence require increased rates of flow to cause filling. 

The above comments on the filling mechanisms are the results of the 
writer’s observations on a barrel length substantially longer and with an 
appreciably greater barrel diameter than any used by the authors. 

Owing in part to the inferred importance of air entrainment phenomena on 

“some of the filling mechanisms and in part to previous experimental work 
which demonstrated the substantial effect of approach flow turbulence level 
on the filling mechanism of rounded and pipe socket inlets at lengths of 14 to 
41 diameters, additional tests were made with the square edge inlet with a 
63.9 diameter long barrel in which the approach flow turbulence level was 
varied by means of screens and grids. The width of the approach channel 
was not varied, being 13.1 diameters in all cases and the grid and screens 
were installed 6 diameters above the inlet and 10.4 diameters downstream of 
the gravel baffle used in all the tests. 

These data are shown in figure 13. The points plotted are the averages of 
‘as many as 5 consecutive runs under the same conditions. It will be noted 
in figure 13 that two series of tests were made with screens. In the first 


A. G. Anderson and C. E. Bowers, Project Report No. 37, St. Anthony 
Falls Hydraulic Laboratory, University of Minnesota. Appendix V, 
Table II, runs 460 and 461. 
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series of tests one double screen was used. This consisted of a 1-inch by 
4-inch wooden frame on each side of which commercial insect screen was 
tacked. This screen assembly thus consisted of two screens placed 4 inches 
apart. In the second series of screen tests two double screen assemblies 
were used resulting in 4 screens placed 4 inches apart. The insect screens 
used were made up of wires of 0.011-inch diameter spaced 14 wires to the 
inch in one direction and 16 to the inch in the other direction. It will be ob- 
served that there is no significant difference in results between the two 
series of screen tests, from which it is inferred that the screens used re- 
duced the approach flow turbulence level sufficiently so that any further re- 
duction in level would not significantly affect the filling mechanism. 

The grid used consisted of vertical 2-inch wooden slats placed 1 inch 
apart. The fact that the data for the grid lie between those for the screens 
and the gravel baffle would indicate as noted previously by Rouse? that the 
flow from a crushed stone baffle of high solidity ratio may of itself introduce 
a high level of turbulence. In the case of the writer’s channel, the turbulence 
level owing to the gravel baffle was apparently sufficiently high so that the 
effect of the grid, of the size used was to reduce rather than increase the 
turbulence level at the inlet. 

It will be noted from figure 13 that the filling mechanism over the range of 
slopes above 0.6% is sensitive to approach flow turbulence level and that a 
substantial effect exists. 

Of significance also is the fact that type B filling described above was not 
observed in the tests with screens, and it would appear that approach flow 
turbulence level affects regime change from sluice type flow to full conduit 
flow in at least two ways. First, it may exert its effect of producing full con- 
duit flow at a lower H/D simply by accentuating the intensity of a filling 
mechanism, for example, the air entrainment process at the greater slopes, 
or it can exert its effect by completely altering the type of filling mechanism 
occurring at a given Slope, as was the case for type B filling with screens in 
the approach channel. 

The data described above on the five distinct filling phenomena observed it 
tests with a 63.9 diameter long culvert are believed directly pertinent to the 
authors’ results for shorter conduits only in regard to the general applicabili 
ty of equation (14). 

The authors’ initial relation with reference to the development of equation 
(14) is equation (5) which is stated to follow from dimensional considerations, 
This is true, of course, only if it is assumed that at filling 


f (L, D, g, Q, n, S) = 0. (25) 


In view of the demonstrated effect of approach channel turbulence, the im- 
portant part played in the filling mechanism by vortex action and air entrain- 
ment in type B filling, and the existence for the higher slopes of a filling 
phenomenon involving a substantial air pressure gradient accompanied by 
violent air entrainment and strong local circulation of air in the inlet region, 
as observed in the writer’s tests with the longer barrel and larger barrel 
diameter, the assumption of general applicability of equation (25), asa 


9. An Investigation of Flow Through Screens, by W. D. Baines and E. G. 
Peterson, Transactions of the American Society of Mechanical Engineers, 
July 1951, discussion by Hunter Rouse. 
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| aie to the derivation of equation (14) would appear to be an assump- 
tion of importance which deserves experimental exploration. 

The authors reduce by one, the number of dimensionless groups of varia- 
bles in equation (5) or (7) by use of the gradually varied flow equation and 
arrive at equation (14). This procedure involves the further assumption that 


2. : Aree | 

since the variables L/D, n'/D /3 and S may be combined to the groupings 

1/3 

ob SD 

SS and 5 to represent backwater effects upon the filling mechanism 
n 


that a Similar combination of L/D, pitt iat and S to represent (referring 
again to the writer’s observations on the filling phenomena in longer and 
larger barrels than used by the authors) the two fluid phenomena of air- 
carrying vortices and the complex air entrainment process of types B, 2, and 
3 filling is also valid. That is, there appear to be three basic physical phe- 
nomena controlling the five distinct filling mechanisms over the range of 
slopes tested by the writer. These involve (a) strictly backwater effects, 

(b) vortex action over the inlet supplying air to the culvert barrel, and (c) a 
violent air entrainment process near the inlet accompanied by a substantial 
adverse air pressure gradient over the free surface, and an energy dissipa- 
tion process involving the transport of air from a region of low pressure toa 
region of high pressure. 

The authors in deriving equation (14) from equation (5), on the basis of 
tests with barrels appreciably shorter and smaller in diameter than the 
writer’s, base their development of equation (14) on phenomenon (a) only. 
The validity of this development with reference to phenomena (b) and (c) 
would appear to constitute an important assumption requiring some 
justification. 

The data of figure 13would simply, regardless of the scaling characteris- 
tics of phenomena (b) and (c) that first, equation (14) requires an additional 
term representative of approach flow turbulence characteristics, and second 
that the data of figure 7 is applicable only to installations having approach 
flow turbulence characteristics similar to the authors’ experimental channel. 

The difficulty in establishing a model law for the square edge inlet without 
full consideration and control of turbulence level in the approach flow is 
illustrated by the data of figure 14, in which a portion of the authors’ data of 
figure 7, together with the writer’s data of figure 13 are plotted in accord- 
ance with the relation 


1/3 
H L 2 BDLey 26 


which follows directly, of course, from the authors’ equation (14). In the 


2 4, 
case of the writer’s data the value of pi/ 3m used was 1.24 x 10 which was 
computed from a full flow, Darcy-Weisbach resistance coefficient determina- 
tion of 0.0149 for the zone of established flow at a Reynolds number of 2.90 x 


10°, which in the writer’s tests corresponded to a value of @/p! 2 of 8. If it 
were assumed that the difference in test results between the two sets of data 
in figure 14 were due only to difference in turbulence level between the two 

channels it would follow that the turbulence level in the authors’ channel was 
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lower than in the writer’s. The difficulty in accepting this interpretation of 
the data of figure 14 stems from two facts. First, as noted in figure 13 the 
writer’s tests with screens in the approach channel consisted of two series c 
data. In the first series, one double insect wire screen was used, and in the 
second series, two double insect wire screens were used. As indicated in 
figure 13 no consistent difference in results was noted between the two serie 
of data and it is to be inferred that the turbulence level in the writer’s tests 
with two double sets of screens was sufficiently low so that any further re- 
duction in approach flow turbulence level would have no effect upon the fillin 
mechanism, and hence that the difference in test results is not attributable 
to the possibility of a lower turbulence level in the authors’ channel than in 
the writer’s. 

This conclusion would appear to be further reinforced by the fact that the 
authors’ gravel stilling baffle was located (as determined by scaling from 
figure 1) approximately 8 barrel diameters from the inlet while in the writ- 
er’s setup the gravel baffle was located 16.4 diameters from the inlet, which 
from consideration of the work of Baines and Peterson! with grids of high 
solidity ratio, would lead to the conclusion that perhaps the turbulence level 
in the authors’ channel was higher than in the writer’s, rather than the 
reverse. 

However, the dangers of drawing general conclusions from the apparent 
conflict of the two sets of data in figure 14, plotted in accordance with equa- 
tion (26) (or equation 14) are obvious. First, the authors’ tests are for a 
non-supported efflux jet, while in the writer’s tests an apron support was 
used. Secondly, the writer’s data were obtained with a headwall 1.33D high 
with modeled spill cones and with a modeled embankment slope, while in the 
authors’ tests, referring to the data of figure 7, spill cones were not used 
and the headwall was the full depth of the water in the upstream pool, withou 
any modeled embankment slope. Third, the method used by the authors in 


tg AiO 
determining n /D /3 is unknown. If the authors’ determination was compute 


on the basis of including the pressure drop in the zone of flow establishment 
their computed value will not, of course, be directly comparable to the 
writer’s value for the zone of established flow. 

It can only be concluded that further consideration of the problem is 
desirable. 

The writer’s comments have stemmed principally from the experimental 
data given in figure 13 on the effect of turbulence level on the filling or prim 
ing phenomena in a culvert barrel, larger in diameter and appreciably longe 
than any used by the authors. The purpose of these comments has been to 
further illustrate the complexity of a flow problem, which from its boundary 
geometry could be assumed to be simple. 


10. An Investigation of Flow Through Screens, by W. D. Baines and E. G. 3 


Peterson, Transactions of the American Society of Civil Engineers, ~ 
July 1951. 
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Discussion of 
“PROPORTIONAL WEIRS FOR SEDIMENTATION TANKS” 


by J. C. Stevens 
(Proc. Paper 1015) 


GUIDO DI RICCO.!—The problem considered by the author is part of the 
more general problem of finding the profile equation x = f(y) for a weir, in 
which the discharge equation Q = F(h) is given. This general problem is ex- 
yressed by the equation 


h 
eave Neate ae (1) 
° 


n which c is the discharge coefficient, either previously known or unknown, 
ind is either a constant or a function of h. 

Such a problem has been solved2 and the following general formula was 
ybtained: 


ie (2) 


ES), 


there were also shown? the particular derivations of known equations such as 
he Pratt formula for a linear or proportional weir (a formula established by 
ratt in a different and rather complicated way, i.e., with a series expansion), 
he exponential-equation weir, and others. 

In this manner, for example, for weirs with a given equation, 


Q=c gah m>1/2 3 0 , @ = const. (3) 
ne following expression was obtained from (2): 


a 2m4 1) 0 (m) , om - 3/2 (4) 
Yr (m+ 1/2) * 


rom which one obtains, for m = 1, the proportional weir of Rother. 

Moreover there also was derived from (2) the profile equation for a weir 
aving a fixed percentage error for a fixed head; the “exponential” weir of 
taus and Standen was thus derived. 


, Prof. of Hydraulics, Univ. of Rome, Italy. 
. G. Di Ricco—Equazione di forma di uno stramazzo dedotta dalla sua 
equazione di portata—L’Energia Elettrica, ottobre 1956. 
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It was also shown? that one could determine the profile of a weir with a 
tiven discharge equation, where the dependence of the discharge coefficient 
- on h head is unknown, by (2) and experimental tests. Practical examples 
vere then solved. It appears as though such studies are unknown to the 
uthor inasmuch as no mention of them is made. 

However, much more work does exist. Weir profiles with a polygonal 
‘-ontour which were easier to construct than those with a curvilinear contour 
vere also considered. These weirs were analytically treated as the sum or 
lifference of either triangular or rectangular notches, and the research was 
ipplied for the design of a semi-module for irrigation-flow measurements, 
similar to the cone module, but with a linear discharge equation. 

The weir under consideration was placed at the end of a rectangular flume 
vith a horizontal bottom of width L and shape as shown in Fig. 1; its dis- 
charge equation for head ranging from 2a to 10a and for + 2% precision was 
‘xpressed by 


Q = 2LVa(h + 2 a) (5)5 


Q in wm.per sec; L, a, and hin meters). This weir was also examined 
vyhile submerged and its discharge equation which was a hyperbolic type® 
vas given. : 
The polygonal contour of weirs was applied to the design of a hydrometric 
levice for flow measurements with linear discharge equation Q = h - 1; such 
. device is now being used at the Istituto di Idraulica dell ’Universita di 
toma” (Fig. 3). 
Recently, the general problem considered herein has been treated in an 
asier way by the use of polygonal-profile weirs which could fit to any form 
f the discharge equation, given either analytically or graphically, upon heads 
n the arithmetic progression a, 2a, 3a, . . ; the value a represents the lower 
imit of validity of the equation. 
That is, given the value of Q for a, 2a,... , na, and letting Q, = Q(ra), the 
rofile of the weir is defined by the following width L, at the variable level 
yr =a: 


r=1 
L,. = eee b, ae + at, ’ (6) 


n which 


. G. Di Ricco—Stramazzi con data equazione di portata—L’Energia Elettrica, 
febbraio 1939. Lal 

G. Di Ricco—Edificio di misura ad equazione di portata lineare—Rivista 
del Catasto, 1940, n' 3 e 4. 

For L = 0,5 meters and a = 0,04 m., (5) becomes Q = 0,2 h + 0,005. 

G. Di Ricco—Edificio di misura ad equazione di portata lineare. Funziona- 
mento sotto rigurgito—Rivista del Catasto, 1949, n°1. 
B. Gaddini—Costruzione di uno stramazzo ad equazione lineare di portata— 
-Rivista del Catasto, 1951, n°3. ri 

G. Di Ricco—Stramazzi a contorno poligonale—Rivista del Catasto, 1954, 
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Table I 


0 150,8654 5732270 
1 = 1,5000 = 313,5023 - 118,9206 
2 350711 65154728 247,1334 
3 - 6,3771 - 1353,7864 = 51959292 
4 13,2524 10671408 
5 - 72,5965 — 27,5386 


If cy = c(ra) is either a known constant or a known function, the problem is 
immediately solved; if such is not the case, the profile must be found by a 
“trial and error” method, first assuming approximate values of Cy and test- 
ing these on the corresponding weir. 

In (6) the value Lo of the base width is arbitrary, such freedom is useful 
to give the weir an aesthetic shape. 

Two examples of weirs, designed by use of the foregoing, are reported on 
herein. 

The first concerns a cylindrical channel with a given trapezoidal cross- 
section and bottom slope where uniform flow is wanted, whatever the dis- 
charge may be. 

The discharge equation was taken in the Bazin form (metrical system); 


The cross-section is shown in Fig. 4, where the discharge equation and the 
weir profile, too, are indicated. 

The channel was 6 meters wide at the bottom, the slope of sides was 2: 3. 
rhe average loss of head was S = 0,0003, c was assumed to be 0,6, and 
1 = 0,5 meters. 

The second example considers a weir suggested for model studies of the 
Tiber River. The discharge equation Q(h) was numerically known (data of 
‘Ripetta” Hydrometer in Roma), and it was assumed that c = 0,6 a = 2 meters. 

Figs. 4 and 5 show the contours (to a distorted scale) of the two examples 
summarized herein. It should be noted that all measurements are in meters. 

The computation of the last shown method is not difficult or impossible, 
‘egardless of the given equation, either Q(h) or V(h), and is resolved in a few 
ninutes only by the use of Table I. This allows one to work with contours in 
umber of 10 linear reaches for each side; usually, however, 4 or 5 reaches 
hall be sufficient, due to the allowance that the weir might not reproduce its 
quation at lowest heads. 
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Discussion of 
“THE APPLICATION OF SEDIMENT-TRANSPORT 
MECHANICS TO STABLE-CHANNEL DESIGN” 


by Emmett M. Laursen 
(Proc. Paper 1034) 


SAM SHULITs,! M. ASCE.—Mr. Laursen summarizes compactly, care- 
fully, completely and in an illuminating manner just how a hydraulic engineer 
should proceed to design a stable channel with the knowledge of today. At the 
same time his cogent comparison of available methods and formulas shows 
clearly and precisely the difficulties and uncertainities with which stable- 
channel design is beset. 

In the last quarter-century, circa 1930 to the present, research in 
sediment-transportation has been increasingly active and inspired. Particu- 
larly true has this been in the United States in the last 10 to 15 years. It is 
then startling to be reminded by the author, first, that substitution of the 
velocity at the edge of the laminar subiayer (vg) for the velocity at the level 
of the moving grain (vg) is tantamount to the 1879 deduction of P. duBoys, 
that q. is a function of T; and second, that transformation of the duBoys bed- 
load formula yields the 1926 permissible canal velocities of Fortier and 
Scobey. To many hydraulic engineers, it may be a strained application to ap- 
ply at the open-channel bed the laminar sub-layer velocity equation based on 
Smooth pipe experiments, whether or not the laminar sub-layer is actually 
existent. The only conclusion is that research in open-channel hydraulics 
and in sediment-transport still has not found, perhaps has not sought enough, 
the appropriate fundamental concepts of fluid mechanics. Until more particu- 
larized facts for open channels are available, one can only use, in the analyti- 
cal process of fluid mechanics, what is on hand from the fields of aerody- 
namics and pipe turbulence. Unfortunately, the results are occasionally the 
old tools with new handles, at least in the two instances just cited. 

Comparative studies of bed-load formulas have been made at the Hydrau- 
lics Laboratory of The Pennsylvania State University. One result from this 
work (reference A)* is presented in Table A, to illustrate the author’s state- 
ment that the difficulties in stable-channel design lie in “assigning numerical 
values to coefficients and even in choosing a formal equation.” In addition to 
the tabulated formulas, there exist today at least two more, the Meyer-Peter 
and Mueller formula (B) and a modified Einstein procedure (C,D). This 
jungle of formulas is so confusing that the engineer confronted with a sedi- 
ment-transportation problem usually resorts to the choice of one or two of 
he most widely quoted in the publications crossing his desk, accompanied by 
snawing uncertainties. 

In the table X and Y are coefficients varying from formula to formula, 
ind sometimes within one formula; T is the tractive force and Tp the 


. Associate Prof., Dept. of Civ. Eng., Pennslyvania State Univ., University 
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inceptive or critical tractive force; n the Manning or Ganguillet and Kutter 
roughness coefficient; V the mean velocity and V, the critical mean velocity; 
m an exponent varying from formula to formula; q the discharge and qo the 
critical discharge, both per unit width; S the slope; d the grain diameter, R 
the hydraulic radius; Y p and YF the densities of bed material and fluid 
(water), respectively while f is the mathematical symbol for “a function of.” 

Though Table A does incorporate an attempt to demonstrate similarities, 
these resemblances in no way imply even a similar range of quantitative re- 
sults. Some of the formulas contain unsuspected sensitivities which vitiate 
simplifying assumptions. For example, a cursory glance at the Schoklitsch 
formula for the inceptive tractive force for sand (E), 


might lead to the conjecture that the 7/6 power of S is an exaggerated refine 
ment—that is, that the lst power might be an adequate approximation. Yet 
there is an appreciable difference between S7/6 and S: at S = 0.00001, S$7/6 = 
0.1468 and at S = 0.01, S7/6 = 0.4648. 

Even for those bed-load formulas where the basic data, usually from 
flume tests, cover the same range, the computed results may be one to two 
hundred percent or more apart. These remarks pertain particularly to the 
DuBoys-Straub, Waterways Experiment Station, Schoklitsch, Casey, Haywooc 
and Meyer-Peter formulas in the comparisons made at the Pennslyvania 
State University. Though the occurrence of such disagreements is well- 
known, the array of formulas in the author’s Tables I and II and the writer’s 
Table A underlines the difficulties of the formula dilemma. 
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} ‘$ 1019(ST4), 1020(ST4), 1021(ST4), 1022(ST4), 1023(ST4),1024(ST4)°, 1025(SM3), 1026 
M3), 1027(SM3), 1028(SM3)¢, 1029(EM3), 1030(EM3), 1031(EM3), 1032(BM3), 1033(EM3)°. 


GUS f: 1034(HY4), 1035(HY4), 1036(HY4), 1037(HY4), 1038(HY4), 1039(HY4), 1040(HY4), 
D41(HY4)©, 1042(PO4), 1043(P04), 1044(PO4), 1045(PO4), 1046(PO4)°, 1047(SA4), 1048 


4)€, 1049(SA4), 1050(SA4), 1051(SA4), 1052(HY4), 1053(SA4). 


BER: 1054(ST5), 1055(ST5}, 1056(ST5), 1057(ST5), 1058(ST5), 1059(WW4), 1060(Ww4), 
ne seatyTid) 1063(WW4), 1064(SU2), 1065(SU2), 1066(SU2)°, Toasters)’: 1068 
, 1069(Ww4). 


E OBER: 1070(EM4), si ae aan ey aaa ed eens pees, 

WOrT(HY5), 1078(SAd), 1079(SM4), 10B0(Sm4), 1081 (SM), 1082 ; 
rae 1087(SA5), 1088(SA5), 1089(SA5), 1090(HW3), tos1(EM4)°, 1092 
, 1093(HW3)°, 1094(PO5)°, 1095(SM4)°. 

§ [BEF eae 109 1098(ST6), 1099(ST6), 1100(ST6), NOM(ST6), U02(IR3), 1103 
ne Leino, tee 6(ST6), Coreen, 1108(STS), 109(AT3), 110(AT3)°, uudR3)*, 


WR: 1113(HY6), (HY6) 5 SA6), 1117(SU3) 1118(SU3) 1119(WW5), 
v3 a fies ww) ey erat Sa areca 1125(BD1)°, 1126(SA6), 1127 


: Pio a Teas 1132(P08), 1133(PO8), f194(P08), 
VOLUME 83 (1957) 
; 1137(CP1), 1138(2M1), 1 eveneey 1140 


i. asa) oe ee! 


, 1141(8M1), 1142(SM1), 
1), 1149(ST1), 11 yy 86 
seen; SEEN fs 


i Hoard, nels usa, 


ees: 


- 


es 2. 
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